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I. INTRODUCTION 



A. PURPOSE 

This study continues the work started by Major D. D. 
Tucker [Ref. 1] and continued by Major John R. Robinson 
[Ref. 2] in their respective theses submitted at the Naval 
Postgraduate School (NPS) in September 1985 and March 1986. 
Their works [Refs. 1,2] dealt with obtaining better attri- 
tion rate estimates for the Marine Corps officer manpower 
model than the ones currently in use. Both Tucker and 
Robinson in the "Recommendations" section of their theses 
[Ref. 1: p. 72] [Ref. 2: p. 69] stated that further work on 

how to update the estimates from year to year was needed. 
This study investigates this "yearly updating" problem. 

The primary purpose of this work is to investigate the 
efficacy of the exponential smoothing model as a yearly 
updating model for Marine Corps attrition rates. The reasons 
for choosing the exponential smoothing model for study are 
outlined in Chapter II. 

A secondary objective of this study will be to compare 
the performance of the attrition rate estimators introduced 
by Robinson [Ref. 2], which are introduced later in this 
chapter, to the current Marine Corps' estimator. 

B. BACKGROUND 

Much of the background detail for this study is well- 
documented in previous works. Tucker [Ref. 1: pp. 15-38, 

128-138] explains the Marine Corps Officer attrition and 
promotion structure, the structure of the data base used by 
the Marine Corps in its officer planning process, the 
process itself, and how the data from the Navy Personnel 
Research and Development Center was transferred to the NPS 
computer system. Robinson's study [Ref. 2: pp. 11,14-20] 
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contains a complete summary of the estimation methods and 
aggregation procedure he used. For continuity, these estima- 
tion methods and this aggregation scheme are used in the 
present work. 

Robinson's thesis [Ref. 2: pp. 74-75] also includes an 
explanation of the Freeman-Tukey arcsine transformation, 
which is used to stabilize the variance of the empirical 
loss rate estimates. The empirical process is assumed to 
have a binomial( n,p) distribution with parameters 

n = central inventory for the year 

p = probability that an inventory unit leaves the 
Marine Corps during the year. 

The Freeman-Tukey arcsine transformation provides a second 
scale for comparison of the estimators within which the 
estimates have a more stable variance, and additionally, 
behave more like a normal distribution (see Appendix B). 

Both Tucker and Robinson [Refs. 1,2] used data from 
years 1977-1980 to obtain their loss rate estimates, and 
then used years 1981-1982 to validate them. Robinson 
[Ref. 2] also had year 1983 to use for validation of his 
estimates. Their estimates [Refs. 1,2] worked fairly well in 
predicting attrition one year into the future, but their 
success faded noticeably as the they tried to predict two or 
three years into the future. Thus there appears to be a 
time-varying component in the attrition rates. This leads 
to the problem of updating the loss rate estimates each year 
in order to better forecast future attrition. 

C. PREVIOUS WORK 

Major Tucker [Ref. 1] showed that the James-Stein 
shrinkage estimator was better than both the current method 
used by the Marine Corps and maximum likelihood estimators. 
He also showed that the James-Stein technique will provide 
estimates for those small cells which have no attrition. 



11 



i. e. , those cells whose maximum likelihood estimator must 
be zero [Ref. 1]. Also, Tucker stated in his summary that 
the small inventory cells present a problem in loss rate 
estimation because "some of the cells are empty for struc- 
tural reasons while others are empty by chance" [ Ref. 1: p. 
70] . 

Major Robinson [Ref. 2] tried to combat the "small cell" 
problem, i. e. , the problem with estimating attrition rates 
for those cells with low inventory figures, by introducing 
the limited translation James-Stein technique of Efron and 
Morris [Refs. 3,4]. He showed that this technique improves 
upon the James-Stein estimates used by Tucker [Ref. 1] in 
estimating the rates for small cells [ Ref. 2] . Robinson 
[ Ref. 2] also introduced a transformed scale cell average 
(TSCA), an estimator corresponding to zero shrinkage in the 
James-Stein technique, which in many cases outperformed all 
other estimators. 

D. AGGREGATION METHOD 

The United States Marine Corps Officer Corps numbers 
approximately 20,000. Each officer below the rank of briga- 
dier general is cross-classified into one of 40 military 
occupational fields (OF), 31 lengths of service years (LOS), 
corresponding to 0 to 30 years in the Marine Corps, and 10 
grades, from warrant officer 1 to colonel, for a total of 
12,400 categories. 

Almost half (6149) of these categories, or cells, are 
"structural zeroes" in inventory. "Structural zeroes" occur 
due to Marine Corps policy concerning promotions and because 
certain combinations of OF, LOS, and grade never occur, 
e. g. , there are no colonels with just 2 years in service. 
These cells exist only in theory, not in practice, and are 
therefore not included in the feasible region of cells. 

A vast majority of the 6251 feasible cells are low 
inventory cells. Because of this, it is quite difficult to 



12 



obtain useful stable attrition rate estimates for these 
cells and it is wasteful to try to treat each cell individu- 
ally. There is much communality of behavior among clusters 
of cells, and the grouping of cells into aggregates of like 
characteristics can ease the bookkeeping burden as well as 
provide the desired stability. Ideally, aggregation schemes 
can be found for which the aggregates behave in a predic- 
table manner and for which meaningful conclusions may be 
drawn from statistical tests. However, current Marine Corps 
practice groups the cells according to organizational and 
operational considerations, producing aggregates that will 
not necessarily conform to any specific statistical 
behavior. 

H. Amin Elseramegy used the CART algorithm to find 
aggregations with predictable statistical behavior with 
encouraging results in his thesis submitted at NPS in 
December, 1985 [Ref. 5]. His results cannot be regarded as 
definitive because of operational considerations ( e. g. , 
excessive computer running time required some pre- 
aggregation), but can serve to guide future work. 

The current Marine Corps manpower model places all occu- 
pational fields into four categories: aviation (OF 72, 75), 

combat support (OF 13, 25, 35), ground combat (OF 03, 08, 

18) and other (includes 32 occupational fields). This 
aggregation scheme is used by both Tucker and Robinson 

[Refs. 1,2] and will be used again in the present work. For 

continuity, the aviation category will include only OF 75 as 
it did in both of their works [Refs. 1,2]. 

E. ESTIMATION METHODS 

Robinson in his thesis [Ref. 2] compared six loss rate 
estimators. These estimators are: 

1) Original Aggregate ( AGG ORIG) -- the current Marine 
Corps estimation method. The occupational fields are 
placed into the four categories mentioned in the 

preceding section. Past attrition rates from 1977 to 
the present are subjectively weighted and averaged for 
each aggregate. The grand mean of the weighted 
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attrition rates serves as the loss rate estimate for 
all cells (OF, LOS, grade) within the aggregate. Both 
Tucker and Robinson ["Refs. 1,2] found methods superior 
to this one. 

2) Transformed Aggregate ( AGG TRANS) -- computed by 
transforming the empirical attrition rates using the 
Freeman-Tukey equation and then calculating the mean 
of the transformed values within each aggregate. This 
is followed by an inversion to the original scale. 
Again, this is a single number used for all cells 
within the aggregate. 

3) Maximum Likelihood Estimator (MLE) -- Is calculated by 
summing all leavers (over time) in a cell and dividing 
by the total cell inventory (over time) for the 
estimation period. This estimator is the MLE for the 
binomial distribution described in the previous 
section. There are problems with using the MLE in this 
setting. While the estimate is indeed unbiased, it is 
unstable due to the abundance of small cells causing 
the possibility of a wide range of values. Also, this 
MLE assumes that each year represents an identical 

E opulation, while the data shows that a cell's 
ehavior can change drastically within a few years' 
time. A complete discussion of the problems with using 
the MLE can be found in [Ref. 2: pp. 17-18]. 

4) Transformed Scale Cell Average (TSCA) -- computed by 

transforming the cell inventory and loss data ana 

calculating the time average for each cell over the 
estimation period. Inversion of the results provides 
attrition estimates in the original scale. Use of the 
Freeman-Tukey transform holds down the variability of 
this estimator, which was mentioned above as a 
shortcoming of MLE. This estimate performed 
surprisingly well in Robinson s analysis [Ref. 2]. 

5) James-Stein Estimator (JS) -- operates from the basic 

notion that by shrinking estimates toward the grand 

mean, the size of the sum of squared residual errors 
will be lessened. The shrinking is applied to the TSCA 
estimator. An optimal shrinkage factor is found for 
each aggregate, and the cell means are shrunk toward 
the grand mean by that amount. See Figure 1. 1. The 
optimal shrinkage factors used here are those found by 
Robinson in his thesis [Ref. 2: pp. 34-36]. Notice 

that since the shrinkage is done in transformed space, 
the assumptions of normality and stable variance 
required by James-Stein are less compromised. 

6) Limited Translation James-Stein (LTJS) -- intuitively, 
it does not seem quite right to shrink all of the 
values towards the grand mean. After all, extreme 
values do occur occasionally, and their effects should 
be represented in an analysis of the data. This 
estimator deals with this problem by limiting the 
translation of extreme values toward the grand mean. 
From Figure 1.2 one sees that there is an interval 
about the grand mean within which full James-Stein 
shrinkage occurs, while outside this interval, the 
shrinkage is diminished. There is a factor, d, which 
controls the width of this interval. Robinson [Ref. 2: 
p. 37] found the optimal values for this factor, and 
they are used in the present work. For a detailed 
explanation of this estimator, see Robinson s thesis 
[Ref. 2: pp22-25] . 
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F. RESULTS 



Exponential smoothing provides a valid yearly updating 
model for the Marine Corps attrition rates in all cases. 
The estimates produced by exponential smoothing are, more 
often than not, better than those produced by the methods of 
Robinson [Ref. 2] , and require far less data. This model 
shows its extreme flexibility by pointing out an external 
change in the aviation environment that occurred in 1981 and 
actually outperforms the Robinson estimates for this aggre- 
gate. 

In this study, as in Robinson's [Ref. 2], the TSCA, JS, 
and LTJS estimators forecast attrition rates better than the 
current method used by the Marine Corps. The maximum likeli- 
hood estimator performs well in the transformed scale, but 
has extremely large errors in original scale in some cases. 
This is due, in part, to the unsuitability in these cases of 
using the optimal smoothing constant values for transformed 
scale to produce smoothed estimates in original scale for 
this estimate (see Appendix C). However, as in Robinson's 
work [Ref. 2], no estimation method emerges as the clear-cut 
"best" choice. 

The ability of exponential smoothing to update attrition 
rates is encouraging. However, having studied such a limited 
sample of data, we believe that its implementation should be 
delayed until further studies can be done when more data 
becomes available. 
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II. YEARLY UPDATING METHODS 



A. GENERAL 

A yearly updating model is one which allows newly 
received data to be combined with data from the past to 
update the estimate produced by that past data. By updating 
the estimate, forecasts produced by the prior estimate 
should be improved upon. 

B. CRITERIA FOR MODEL SELECTION 

There are many methods available to update parameter 
estimates, such as attrition rates, as more data becomes 
available. To choose among these methods, Robert Goodell 
Brown [Ref. 6] suggests using the following criteria: accu- 

racy, simplicity of computation, and flexibility of rate of 
response. Of the three, only simplicity of computation is 
not a major concern for the yearly updating model. In 1963, 
when Brown published his book, simplicity of computation was 
important for a model because of the relative inefficiency 
of computers of that time as compared with those of today. 
Since today's computers are millions of times faster, we 
will not be concerned with this criterion in our model 
selection. However, a criterion which is important in 
choosing a yearly updating model is the amount of data that 
needs to be stored in order to produce the estimates. 

Naturally there are other criteria which may be used to 
select a parameter estimation model. However, the three that 
seem most appropriate in determining the yearly updating 
model used in this study are: 

1) accuracy 

2) size of data 

3) flexibility of rate of response . 
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1. Accuracy 

Fulfilling the Marine Corps' need for more accurate 
attrition rate estimates is one of the primary purposes of 
this pilot study. More accurate forecasts are produced by 
such estimates, thereby avoiding some of the costly overages 
and underages in inventory which result from the use of the 
current estimation system. 

One of the purposes of finding a yearly updating 
model is to improve upon the forecasts of Tucker and 
Robinson [Refs. 1,2], for which accuracy dropped off sharply 
after the first year. The yearly updating model selected 
should be able to forecast attrition rates one year in 
advance at least as well as Tucker and Robinson's 
[Refs. 1,2] did, and also improve greatly upon forecasting 
rates two or three years from the present. 

In reality, it is the estimation of attrition rates 
two years from the time of the most recent data that should 
be of major concern. Loss and inventory data for a fiscal 
year is not generally available until halfway through the 
following fiscal year. Therefore, with the data available, 
the next time period in which attrition rates need to be 
forecasted is the year following the year in which the data 
is received, which is two years ahead of the most current 
available data. Thus the accuracy of the estimates for two 
years from the time of the most current data, for which 
Tucker and Robinson [Refs. 1,2] had little success, is one 
of the key figures with which the yearly updating model need 
be concerned. 

2 . Size of Data 

This criterion is very important in choosing a 
yearly updating model. With 6251 feasible cells, it becomes 
unwieldy and costly in terms of computer storage space to 
log year after year of loss and inventory data. Thus, a 
model is sought which can use short summaries of the data 
and still produce valuable estimates. 
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Flexibility to the Rate of Response 

In forecasting, when the current observation is 
different from the forecasted value, there are two possible 
explanations: random fluctuation; change in the pattern of 
the data. If the error is a random fluctuation in the data, 
then the forecasting technique should smooth out the fluctu- 
ation. In order for the model to do this, it should produce 
estimates based on a great deal of past data. However, if 
the error is due to a change in the pattern of the data, 
then past data is rendered irrelevant. The estimate should 
reflect only the recent processes. 

Changes in Marine Corps policies in handling its 
officer corps, which occur from time to time, often cause 
corresponding changes in attrition patterns. Given informa- 
tion about such changes, one should be able to detect a new 
pattern in the loss and/or inventory data resulting from 
them. Thus, we want our yearly updating model to be able to 
smooth out random fluctuations in the data during times when 
the attrition process is stable over a period of years, yet 
still be able to respond rapidly to new conditions. The 
model must be able to easily adjust the number and relative 
value of past observations in producing the estimate, using 
a fairly long series of data for an unchanging process and 
only the most recent observations when a change in the 
process occurs. 

C. SELECTION OF THE YEARLY UPDATING MODEL 

1. Candidate One -- Tucker /Robinson 1 s Methods 

Both Tucker and Robinson [Refs. 1,2] showed that the 
alternate estimators they introduced outperform the aggre- 
gation methods currently employed by the Marine Corps. They 
used four equally weighted years of data, 1977-1980, to 
produce their estimates [Refs. 1,2]. The question that 
arises from their findings is why not use the functions that 
Tucker and Robinson wrote [Refs. 1,2] to compute estimates 
based on five, six, seven, etc. , years of data? 



20 



The answer to this question lies in the criteria for 
model selection outlined in the previous section, of which 
the continuing use of Tucker and Robinson's methods over 
time meet none. The accuracy of Tucker and Robinson's esti- 
mates [Refs. 1,2] dropped off sharply after the first year 
in most cases, and the estimate of attrition two years off 
is very important, as mentioned earlier'. The amount of data 
used to produce their four-year estimates was tremendous; 
storing all of that data plus that of additional years would 
be totally inefficient. Finally, using their methods with 
more data is not responsive to changes in attrition 
patterns. Old data, equally weighted with recent data and 
which may no longer be relevant, is used in this technique 
to produce the loss rate estimates. 

Therefore, using Tucker and Robinson's methods 
[Refs. 1,2] year after year on all of the data available 
since 1977 is not a very good alternative. However, their 
methods are valuable for providing a base estimate for the 
exponential smoothing model to be discussed later. 

2 . Candidate Two -- Moving Average 

A moving average is simply the sum of the most 
recent N observations divided by N. For example, imagine 
that a basketball player scores 18, 15, 25, 22, 20, and 14 
points in his first six games of the season. Let the scores 
be denoted by x^, i = 1,2,..., 6. His scoring average is the 
sum of these six scores divided by six, or 19. 0 points per 
game. Call this Mg. Suppose he scores 24 points in his 
seventh game. Then, if the practice of finding his scoring 
average over the past six games is continued, his six-game 
moving average. My, can be computed by adding the scores 
from the most recent six games and dividing by six. Another 
way would be to subtract 1/6 of the score he achieved six 
games ago, or game one, and adding 1/6 of his total for the 
most recent game, game seven. This produces the six-game 
moving average 
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M 7 = M 6 + ( x 7 ” x 1 )/N 

= 19. 0 + ( 24-18J/6 = 20. 0. 

Using this procedure to find attrition rates would 
be straightforward; one would just find the actual rates for 
the past N years and average them. This method does away 
with some of the problems with data size encountered with 
the previous candidate. Whereas continued use of Tucker and 
Robinson's methods [Refs. 1,2] required the storage of all 
of the inventory and loss data from 1977 forward, this model 
will compute the rates each year for aggregates of cells, 
and only these rates, which number N times the number of 
aggregates, need be stored. This method may also prove to 
be more accurate than the continual use of Tucker and 
Robinson's methods [Refs. 1,2] since data from more than N 
years ago, which may bear little resemblance to current 
data, is eliminated. However, it is difficult to change the 
rate of response using a moving average. If a change in the 
underlying distribution of the data occurs, it will take N 
years for the moving average to fully reflect this change. 
One might suggest to keep N small so that it will respond 
more quickly to changes, but by doing this, the greater 
accuracy produced by larger data sets is sacrificed. Thus, 
N must be chosen so as to compromise between these 
conflicting objectives. 

The moving average is an improvement over the 
previous candidate, but it requires more data and is not as 
flexible as the candidate which follows, the exponential 
smoothing model. However, an understanding of it is helpful, 
as the exponential smoothing model is itself a type of 
moving average. 

3. Candidate Three -- Exponential Smoothing 

It is a disadvantage of the moving average that it 
has to carry all of the rates needed to compute it, albeit 
is a great improvement over the amount needed for the 
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continued use of Tucker and Robinson's methods [Refs. 1,2]. 
Exponential smoothing cuts back even further on the amount 
of data needed. Let us see how by continuing the example of 
the previous subsection. 

Suppose now that after recording the 24 points 
scored by the basketball player in game seven, it is discov- 
ered that all records of the previous six scores have been 
destroyed, but the moving average. Mg = 19.0, still remains. 
If the value of x^ was known. My could be computed. The best 
estimate we have for x^ is that is was equal to the average. 
Mg = 19.0. Using this estimate for x-^, a new estimate of the 
six-game scoring average can be computed: 

Mhaty = Mg + (xy-Mg)/6 

= (l/6)x 7 + ( 5/6 )Mg = 19.833 

Mhaty is an estimate of the moving average My, and 
is called the smoothed value of the average. Mhat^ will 
hereafter be referred to as S^, S standing for "smoothing." 

If the equation used to find Mhaty = Sy is used to 
find each succesive estimate, the definition of the smoothed 
function of the observations is [ Ref. 6: p. 101] : 

S t (x) = a x t + (l-«)S t _ 1 (x) (2.1) 

The smoothing constant, a, is similar, but not exactly equal 
to the fraction 1/N used to find the moving average [ Ref. 6: 
p. 101] . The operation which updates an estimate by adding 
a fraction a of the difference between the current observa- 
tion and the previous estimate to that previous estimate is 
called exponential smoothing [Ref. 6: p. 101]. 

Exponential smoothing discounts past data based upon 
the size of the a parameter. How it does so can be seen by 
substituting for the previous smoothed value the equation 
that produced it from an even earlier smoothed value 
[ Ref. 6: p. 101] : 
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s t (x) = u X t + ( l-«)( a x t _ 1 + ( l-«)S t _ 2 ( X) ) (2.2) 

= a x t + a ( l-«)( x t _ 1 ) + a (l-u)^ x t _ 2 +. . . 

= « Z (l-«) k (x t _ k ) + (l-a)S 0 (x) 

Thus, if the smoothing constant equals .2, then the 
current data point has weight . 2. Previous observations 
have weights . 16, . 128, . 1024, etc. . 

It is seen from the above equations that exponential 
smoothing always requires a prior estimate, to perform 
the update and find S^.. Brown [Ref. 6: p. 102] suggests 
using the simple average of the most recent N observations, 
or M^-l' for the initial value S^.^- In this study, the 
prior estimates are found in a similar manner. The functions 
developed by Tucker and Robinson [Refs. 1,2] and spelled out 
in detail in Robinson's thesis [Ref. 2: pp. 83-110] are used 
to find the estimates over N years of data. An optimal base 
estimate length N will be found for each aggregate using the 
data available. However, so that the results of this study 
may be compared with those of Robinson [Ref. 2], 3-year base 
estimates, corresponding to years 1977-1979, will be used. 
The empirical rates for 1980 are smoothed onto the base 
estimates to produce the updated estimates which can be 
validated on years 1981-1983 just as Robinson's estimates 
were [ Ref. 2] . 

The exponential smoothing model best meets the 
criteria outlined in the previous section. It requires very 
little data to be carried from year to year as compared with 
the other candidates examined; only the last estimates 
obtained by the model need be saved. It is also very flex- 
ible to changes in the pattern of the data. When the 
smoothing constant is small, the function behaves like the 
average of a great deal of past data, whereas large values 
of the smoothing constant allow S(x) to respond quickly to 
changes in the attrition rate process [Ref. 6: p. 102]. Its 
accuracy should be greater than that of the moving average. 
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where the data is equally weighted, since exponential 
smoothing discounts past data. This will allow for more 
recent data to exert greater influence on the estimate, 
which is desirable because the near past generally repre- 
sents the near future better than the distant past. 

Therefore, since it fulfills the criteria better 
than all other candidates examined, exponential smoothing is 
the yearly updating model used in this work. 
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III. EXPONENTIAL SMOOTHING 



A. GENERAL 

Having chosen exponential smoothing for the model to 
update the attrition rates annually, the next step is to 
implement it. The exponential smoothing itself will be 
performed in transformed scale because the variances of the 
estimates are more stable in this scale. Therefore, both the 
base estimates S^_ 2 (x) and observations x^ , which are needed 
to produce the updated estimate S^(x), will be calculated in 
the transformed scale. 

When the loss data is collected by the Navy Personnel 
Research and Development Center, the cell (OF, LOS, grade) 
into which a loss is assigned is the cell to which the 
leaving officer belongs at his time of departure. For 
example, a first lieutenant with 3 years of service at the 
beginning of year i who completes his 4th year of service, 
is promoted to captain, and subsequently leaves the service 
before the beginning of year i + 1 is classified as a loss 
from LOS 4 years in his OF and the grade of CAPT in year i. 
This type of loss classification demands the use of central 
inventory data [Ref. 7: p. 25]. For data grouped by years, 
like the Marine Corps manpower data we have, central inven- 
tory for year i is found by averaging beginning-of-year i 
stocks with beginning-of-year i+1 stocks for each cell. If 
beginning-of-year i+1 stocks are not available, because we 
are at the end of the data set or for whatever reason, the 
central inventory for year i is set equal to the beginning- 
of-year stocks for year i. This occurs in our data for year 
1983, since 1984 is not available for our use. Additionally, 
if losses in year i are greater than the central inventory 
of year i, the central inventory is set equal to the losses 
to avoid the apparent inconsistency of losing more men than 
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you have on hand. Thus, in the above example, the officer 
would have been counted as a 1LT in his OF with 3 years of 
service in the inventory for the beginning of year i. The 
attrition rate computed then, is the central attrition 
rate, which is the number of losses in year i divided by 
central inventory for year i [Ref. 7: p. 25]: 

n»i = Yj/Nj. (3.1) 

The raw data compiled by the Navy Personnel Research and 
Development Center consists of losses and beginning-of-year 
inventories for the 12400 cells mentioned in the 
Introduction chapter. Both Major Tucker and Major Robinson 
[Refs. 1,2] aggregated this data into the four categories 
aviation (OF 75), combat support (OF 13, 25, 35), ground 
combat (OF 03, 08; 18), and other. The central inventory 
matrices were computed and the loss matrices compiled for 
the aviation, combat support, and ground combat aggregates 
for all years available ( 1977-1983 for Robinson and this 
work). The 'other' category was not examined by Tucker nor 
Robinson, nor will it be here. The best summary of the data 
manipulation programs producing the central inventory and 
loss matrices is found in Appendix D of Major Robinson's 
thesis [Ref. 2: pp. 83-104]. In order to be consistent with 
the analyses done in those two works, only the grades of 
first lieutenant ( 1LT) and lieutenant colonel (LTCOL) will 
be examined in this pilot study. Like Tucker and Robinson 
[Refs. 1,2], all 31 possible lengths of service years will 
be included. Additionally, all six estimators explored by 
Robinson [Ref. 2] and discussed in the Introduction chapter 
of this study will undergo exponential smoothing and have 
their forecasting abilities compared. Thus the total number 
of loss rate estimates that will ultimately be computed is 
25284, corresponding to 6 estimators, 7 operational fields, 
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31 lengths of service, and 2 grades. 


The analysis 


will be 


broken down by grade and aggregate, 


as in 


Tucker and 


Robinson's works [Refs. 1,2], thereby 


setting 


up 


6 study 


blocks: 








1) 


Aviation 1LTS 








2) 


Aviation LTCOLS 








3) 


Combat Support 1LTS 








4) 


Combat Support LTCOLS 








5) 


Ground Combat 1LTS 








6) 


Ground Combat LTCOLS . 









The base estimate S^_-^ and empirical estimate x^ are 
found for all six estimates, 31 lengths of service, and 
occupational fields within a study block by using APL func- 
tion ESTIM (see Figure A. 3 ) , developed by Major Robinson 
[Ref. 2: p. 105]. The updated estimate, still in trans- 

formed scale, is found by applying exponential smoothing 
equation 2. 1 presented in the preceding chapter. The 
predicted rates are then validated against actual trans- 
formed figures from the years that follow. How this is done 
is the subject of the next section. 

B. FIGURES OF MERIT 

The figures of merit ( FOM) used by Robinson in his 
thesis [Ref. 2] will also be used here. The basis for the 
decisions concerning the finding of the optimal smoothing 
constants and base period lengths will be the figures of 

merit in transformed scale because the data is so much more 
well-behaved therein. The original scale figures of merit 
will also be calculated and reported in Chapter IV. The 

values of a used to calculate the original scale figures of 
merit will be the optimal values of a found in transformed 
scale. 

We would hope, therefore, that the value of the 

smoothing constant, a, corresponding to the smallest figure 
of merit in transformed scale is very close to that 
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corresponding to the smallest FOM in original scale. 
Fortunately, this is almost always the case. Two exceptions 
are the transformed aggregate estimators for aviation 1LTS 
and ground combat 1LTS, which turn out not to be very 
important due to the overall lackluster performance the 
transformed aggregate estimator turns in throughout the 
analyses. More notably, however, the MLE for combat support 
1LTS, ground combat 1LTS, and ground combat LTCOLS have 
serious discrepancies between values of a producing minimum 
figures of merit in transformed and original scales. These 
differences cast serious doubts upon the correctness of 
using the transformed scale optimal a as the a value in the 
original scale for maximum likelihood estimation. See 
Appendix C for the details of this analysis. 

In transformed scale, the data is approximately distrib- 
uted normally with a stable variance (see Appendix B). Thus, 
a good means of comparison between estimators is a sum of 
squares error (SSE) computation. The sum of squares error is 
defined as: 



SSE = £ (actual-predicted)^ (3.2) 

where "actual" and "predicted" are the actual and predicted 
values for the transformed attrition rate figures for one of 
the validation years (VY). APL function RISKT calculates 
the figures of merit for transformed scale using the above 
equation (see Figure A. 9 ) . 

The original scale does not have a normal distribution 
nor a stable variance for the estimates. The SSE would 
therefore be inappropriate for the FOM calculations in orig- 
inal scale. Robinson [Ref. 2: p. 28] sugggested the use of a 
chi-square statistic for use to compare estimators in orig- 
inal scale. The chi-square statistic is: 
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FOM(t) = y (y it -N it p i ) 2 /N it p i ( 1-Pj^), for all i 



(3.3) 



where and are the losses and central inventory 
counts for the cell in the t*'* 1 validation year, and p^ 
is the inverse transform of the estimate for the i^ cell 
calculated in transformed scale (see Appendix B). APL func- 
tion RISKO (see Figure A. 8), which computes the original 
scale figures of merit, screens out cells with p^ values of 
0 and 1 to prevent the above denominator from having a value 
of zero. Both RISKT and RISKO find the figures of merit for 
all of the validation years available. 

With both of these figures of merit, the smaller the 
FOM, the better the estimates produced. Figures of merit are 
summed over all of the operational fields and lengths of 
service within a study block, thereby producing a single 
number for necessary comparisons. 

C. FINDING THE SMOOTHING CONSTANTS 

As mentioned in the preceding paragraph, it is the 
smallest figure of merit which we seek. Therefore, the 
smoothing constant, a, which produces the smallest FOM for 
an estimator is the optimum « for that estimator. 

In order for the weighting scheme implied by exponential 
smoothing to make sense, a must be between 0 and 1 inclu- 
sive. An estimate with a value of a equal to 0 places no 
weight on the current empirical data point and all of the 
weight on the previous estimate; thus, this value is simply 
the estimate produced by Robinson [Ref. 2]. Conversely, an 
estimate produced with a value of a = 1 is simply the empir- 
ical estimate from the most recent year's data. 

An APL function was developed to produce the figures of 
merit in both transformed and original scale for all valida- 
tion years available and all six estimators for a study 
block for all values of alpha between 0 and 1 by a specified 
step size, which in this study is .02. This function. 
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ALPHAHAT, calculates the base estimates using function 
ESTIM, the empirical data points using function XEOUR ( for a 
three-year base; XFIVE and XSIX for four- and five-year 
bases), which calls ESTIM, and smooths them into the new 
estimates using equation 2. 1. See Appendix A. Arrays of 
figures of merit are produced, which can then be analyzed 
and the optimum a found for each estimator and study block. 

D. LENGTH OF BASE PERIOD 

The length of the base period is important in producing 
the loss rate estimates. If the environment is stable over a 
long period of time, a long base period is preferred to 
smooth out random fluctuations in the data. If the environ- 
ment is in a constant state of turmoil, the base period 
length should be rather small, and frequent updates of the 
base are needed. 

To provide a basis for comparison with the results 
obtained by Robinson [Ref. 2], a three-year base will be 
used in this study, with the fourth year empirical data 
smoothed onto it to produce the attrition rate estimates. 
Thus, years 1977-1980 are used for computation of the esti- 
mates, and years 1981-1983 are used for validation, as in 
Robinson's work [Ref. 2]. 

It should be noted that the choice of three years as the 
length of the base period may not be "optimal," that is, it 
may not provide estimates as accurate as those of a four- or 
five-year base period. Therefore, a short comparison of base 
estimate lengths three, four, and five years is presented 
below. The MOE used is transformed FOM for forecasted rates 
for validation year one. Although we are more concerned with 
the estimators' ability to estimate rates two years in 
advance, this MOE is chosen so that a 5-year base, which can 
only forecast for 1983, may be included. 

To find the figures of merit, the optimal tt, which is 
the value of « producing the smallest FOM for validation 
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year one, is found for each base year length and study 
block. The resulting figures of merit are then compared to 
determine the "best" base period length for an aggregate. 

It should be noted that with the limited number of data 
years we have to work with, the conclusions drawn by this 
short analysis cannot be generalized to future years. The 
length of base period problem should be examined in more 
detail as more data becomes available. 



TABLE 1 

DETERMINATION OF OPTIMAL BASE PERIOD -- COMBAT SUPPORT 



FIRST LIEUTENANTS 



Estimator 


3-Year Base 


4-Year Base 


5-Year Base 




« 


F0M81 


a 


FOM82 


a 


FOM83 


AGG ORIG 


. 50 


1. 863 


1 


1. 669 


. 62 


1. 494 


AGG TRANS 


0 


4. 086 


0 


4. 737 


0 


4. 264 


MLE 


. 54 


2. 402 


. 58 


2. 409 


. 48 


1. 894 


TSCA 


. 52 


1. 610 


. 62 


1. 877 


. 50 


1. 136 


JS 


. 78 


1. 678 


. 80 


1. 817 


. 58 


1. 153 


LTJS 


. 58 


1. 641 


. 66 


1. 861 


. 52 


1. 133 






LIEUTENANT COLONELS 






Estimator 


3-Year Base 


4-Year Base 


5-Year Base 




a 


F0M81 


a 


FOM82 


a 


FOM83 


AGG ORIG 


. 66 


1. 478 


1 


1. 441 


1 


1. 362 


AGG TRANS 


1 


3. 372 


1 


2. 789 


0 


2. 603 


MLE 


. 36 


1. 859 


. 50 


1. 625 


. 40 


1. 638 


TSCA 


. 38 


1. 024 


. 72 


1. 209 


. 60 


. 980 


JS 


. 36 


. 935 


. 76 


1. 124 


. 64 


. 934 


LTJS 


. 36 


. 968 


. 74 


1. 143 


. 62 


. 956 


Note: Figures of 


Merit are 


those 


for transformed 


scale 



1. Combat Support 

Interestingly,' there is an increase in the minimum 
FOM from a three- to a four-year base, with the five-year 
base minimum FOM being smaller than both of them for the 
TSCA, JS, and LTJS estimators for both grades, as well as 
AGG TRANS and MLE 1LTS. The other 4 estimator-grade combina- 
tions show a strictly decreasing trend across all 3 base 
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period lengths. Thus, from Table 1, it appears that a 5-year 
base would be best to use for the combat support aggregate, 
since all transformed figures of merit are smallest for this 
base period length. 



TABLE 2 

DETERMINATION OF OPTIMAL BASE PERIOD -- GROUND COMBAT 



FIRST LIEUTENANTS 



Estimator 


3-Year Base 


4-Year Base 


5-Year Base 




a 


F0M81 


u 


FOM82 


a 


FOM83 


AGG ORIG 


. 58 


2. 976 


1 


3. 066 


1 


3. 933 


AGG TRANS 


0 


19. 842 


0 


20. 923 


0 


24. 811 


MLE 


. 72 


1. 647 


. 82 


2. 127 


. 80 


2. 909 


TSCA 


. 8 


1. 324 


1 


1. 668 


. 90 


2. 559 


JS 


. 92 


1. 460 


1 


1. 664 


. 96 


2. 631 


LTJS 


. 86 


1. 348 


1 


1. 609 


. 94 


2. 573 






LIEUTENANT 


COLONELS 






Estimator 


3-Year Base 


4-Year Base 


5-Year Base 




a 


F0M81 


a 


FOM82 


a 


FOM83 


AGG ORIG 


1 


3. 540 


1 


2. 698 


1 


2. 413 


AGG TRANS 


0 


14. 393 


0 


12. 481 


0 


12. 653 


MLE 


. 42 


1. 389 


. 20 


1. 982 


. 40 


2. 250 


TSCA 


. 58 


1. 510 


. 58 


1. 904 


. 70 


1. 922 


JS 


. 58 


1. 505 


. 60 


1. 832 


. 70 


1. 884 


LTJS 


. 56 


1. 494 


. 58 


1. 881 


. 70 


1. 899 


Note: Figures 


of Merit 


are those for transformed scale 



2. Ground Combat 

From Table 2, an increase is seen from both 3- to 
4-year bases and from 4- to 5-year bases in all cases except 
AGG ORIG and AGG TRANS for LTCOLS. Overall, it appears that 
a 3-year base period produces the smallest figures of merit 
for validation year one, and is therefore "optimal" for this 
aggregate and data set. 

The large values of a for first lieutenants are 
noteworthy. They may represent the inherent variability of 
the process, or they may indicate that the attrition process 
for this aggregate is in an almost continuous state of 
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change from year to year. The reasons are unclear at this 
time, and further analysis with more data is needed to 
determine whether or not this trend persists. 



TABLE 3 

DETERMINATION OF OPTIMAL BASE PERIOD -- AVIATION 



FIRST LIEUTENANTS 



Estimator 


3- 


Year Base 


4-Year Base 


5-Year Base 




a 


FOM81 


a 


FOM82 


a 


FOM83 


AGG ORIG 


1 


3. 976 


1 


4. 253 


. 44 


4. 736 


AGG TRANS 


0 


206. 877 


0 


218. 857 


0 


216. 057 


MLE 


. 8 


2. 177 


1 


4. 287 


. 36 


1. 502 


TSCA 


. 98 


1. 243 


1 


4. 512 


. 38 


1. 176 


JS 


1 


1. 402 


1 


4. 494 


. 38 


1. 277 


LTJS 


1 


1. 386 


1 


4. 518 


. 38 


1. 261 






LIEUTENANT 


COLONELS 






Estimator 


3- 


Year Base 


4-Year Base 


5-Year Base 




a 


F0M81 


a 


FOM82 


a 


FOM83 


AGG ORIG 


1 


7. 898 


1 


10. 959 


1 


6. 992 


AGG TRANS 


1 


24. 733 


0 


31. 473 


0 


27. 999 


MLE 


. 88 


3. 379 


. 66 


6. 637 


. 50 


3. 208 


TSCA 


1 


3. 341 


. 94 


6. 325 


. 7 


3. 698 


JS 


1 


3. 302 


. 96 


6. 311 


. 72 


3. 719 


LTJS 


1 


3. 321 


. 94 


6. 377 


. 7 


3. 684 


Note: Figures 


of Merit 


are those for transformed scale 



3. Aviation 

From Table 3, it is seen that for aviators, the 
4-year base period has a higher FOM than the other two 
candidates in all cases except for AGG ORIG 1LTS, for which 
the FOM for a 4-year base length is in the middle of the 
three values. Thus, the choice for the "optimal" base period 
length centers on the 3- and 5-year candidates. For 1LTS, 
the minimum FOM for the 5-year base is smallest for MLE, 
TSCA, JS, and LTJS, while the 3-year base is "best" for both 
of the aggregate estimators. For LTCOLS, the 5-year base has 
the smallest minimum FOM only for AGG ORIG and MLE, while 
the 3-year base excels for the other four. With six 
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estimators being "best" for both the 3- and 5-year base 
period lengths, an aggregate-wide analysis results in a 
stalemate. Thus, the analysis is broken down by grade, with 
the conclusions being that the "optimal" base period lengths 
are 5 years for aviation 1LTS, and 3 years for LTCOLS. 

The a values in this table, shifting from near 1 for 
a 3- and 4-year base to around . 38 for 1LTS and . 7 for 
LTCOLS for the 5-year base should be noted; they indicate a 
radical change in the data. This change comes about, in 
part, because of the initiation of Aviation Officer 
Continuation Pay (AOCP) in 1981. As explained by Major 
Tucker [Ref. 1: p. 18], AOCP provides a bonus per year to 
aviation officers which in turn obligates continued service. 
The program was applied to all ranks provided the individual 
met certain active duty flight status requirements. This 
action had its desired effect on retaining aviation offi- 
cers, according to the analysis of Major Tucker [Ref. 1: p. 
18] . 

In this analysis, the effect of AOCP is that data 
from the pre-1981 era is not relevant to post-1981 data, 
thereby producing a values of 1. 0 for the 3- and 4-year 
base period lengths for aviation. This means that the esti- 
mates producing the smallest figures of merit are simply the 
prior year's empirical estimates, which makes sense when the 
base years include all pre-1981 data, as they do for the 3- 
and 4-year bases. The values of a in the . 7 range for 3 of 
the 6 estimators for lieutenant colonels and the 5-year base 
indicate that theoretically, 19 r 2 data should be given a 
weight of .7, 1981 data a weight of (.7)(.3) = .21, and 
1977-1980 data a collective weight of (.7)(.3)^ = .09 (see 
Equation 2.2). This is also consistent with the AOCP 
program; the exponential smoothing model is beginning to 
build a base of its own beginning in 1981. This is also seen 
to be happening to the 4 estimators introduced by Robinson 
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[Ref. 2] for first lieutenants; the model is constructing 
its own base. The exponential smoothing model really shows 
its flexibility and overall value in its predictions for the 
aviation aggregate. 

E. STABILITY OF THE SMOOTHING CONSTANTS 

The stability of the smoothing constants is important to 
the validity of the exponential smoothing model. This is 
because if the same value of a produces optimal results in 
the short run, then the exponential smoothing model is a 
good one to use. We say "in the short run" because changes 
in conditions over time that affect the data will demand 
updates in « as well as the base period used. 

Despite the relatively short base periods forced upon us 
by a lack of data, the smoothing constants found for each 
study block appear to be rather stable, i. e. , there are no 
wild fluctuations in the a which produces the minimum trans- 
formed FOM. The few exceptions will be identified in the 
analysis that follows. The stability of the smoothing 
constants will be demonstrated in two ways: within base 
period lengths, and between base period lengths. 

1. Within Base Period Lengths 

The stability of a within a base period length is 
measured by how much a varies in producing minimum figures 
of merit for all of the validation years (VY) available. 
The analysis of the stability of « within base period 
lengths, therefore, is restricted to the 3- and 4-year 
bases, since the 5-year base can only be validated on one 
year, 1983, in our data set. In addition, the two aggregate 
estimators are excluded from this and the following (between 
base period lengths) graphical analyses because of the 
consistently poor performance they show as compared to the 
other four estimators as measured by the figures of merit in 
transformed space seen in Tables 1, 2, and 3. An analysis 
of the stability of a for these two estimators reveals that 
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Figure 3.2 Within 3-Year Base Stability -- Aviation LTCOLS 
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Figure 3.3 Within 3-Year Base Stability -- Cbt. Spt. 1LTS 



39 



/ 





CO CVI 





Figure 3. 4 Within 3-Year Base Stability 



Cbt. Spt. LTCOLS 
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Figure 3. 5 Within 3-Year Base Stability -- Grd. Cbt. 1LTS 
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Figure 3. 6 Within 3-Year Base Stability — Grd. Cbt. LTCOLS 
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Figure 3. 7 Within 4-Year Base Stability -- Aviation 1LTS 
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Figure 3. 8 Within 4-Year Base Stability -- Aviation LTCOLS 
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Figure 3.9 Within 4-Year Base Stability -- Cbt. Spt. 1LTS 
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Figure 3. 10 Within 4-Year Base Stability 



Cbt. Spt. LTCOLS 
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Figure 3. 11 Within 4-Year Base Stability 



Grd. Cbt. 1LTS 
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Figure 3. 12 Within 4-Year Base Stability -- Grd. Cbt. LTCOLS 
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the A GG TRANS estimator shows perfect consistency in its 
optimal a values in all cases, and the A GG ORIG a values are 
also very stable in most cases with no notable exceptions. 

Note that when reading the graphs for the 3-year 
base period length, the solid line is for validation year 1, 
1981, the dot-dashed line is for validation year 2, 1982, 
and the dashed line is for validation- year 3, 1983. The 
graphs for the 4-year base period length have a solid line 
for validation year 1, 1982, and a dot-dashed line for vali- 
dation year 2, 1983. The year corresponding to each type of 
line is seen on one of the graphs, usually the "MLE" graph, 
on each page. 

At first glance, it would appear that these twelve 
graphs show some study cells to have very stable a values 
while others have a values that vary greatly. Before any 
conclusions are drawn, however, attention must be focused on 
the scales of the Y-axes. In almost all cases, these scales 
are so small that they render the curves practically flat 
over the interval covered by the range of optimal a values. 
This indicates stability in a despite the seen difference in 
optimal a values since any a value in the range of optimal a 
values will produce a figure of merit very close to the 
optimal FOM. 

A measure of the disparity caused by the differing a 
values is the maximum percent differential of figures of 
merit produced in the range of optimal a values over all 
validation years. This error is measured for each valida- 
tion year, estimator, and study block by subtracting the 
minimum FOM produced in the range of a from the the maximum 
FOM produced and then dividing by the minimum FOM. The 
maximum of these over the validation years for an estimator 
and study block is the error measure. 

For these figures, this value is generally very 
small, usually between 0 and 2 percent. For a 3-year base 
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period length, however, there are some estimator-study block 
combinations with errors greater than 5% (chosen arbitrarily 
to be an acceptable tolerance), and these are: 

1) aviation 1LTS MLE with 15.84% error 

2) combat support LTCOLS JS with 5.24% error 

3) combat support LTCOLS LTJS with 5.78% error 

4) ground combat 1LTS MLE with 12. 87%' error 

5) ground combat 1LTS TSCA with 8.08% error 

6) ground combat LTCOLS MLE with 5.75% error 

This implies that only 6 of the 24 a values studied for the 
3-year base period may vary too much within validation years 
for that base period length. 

For a 4-year base-period length even better results 
are obtained. The maximum percent differential is the 4.2% 
posted by the James-Stein estimator for combat support 1LTS, 
which also had the largest gap in a values, .24. In fact, 
outside of that study group, only one estimator-study block 
combination, MLE for aviation LTCOLS, has a percent differ- 
ential of more than 1. 1%. Thus, Figures 3. 1 to 3. 12 appar- 
ently show that a is stable within base period lengths for 
both cases we can study using the data set. 

2. Between Base Period Lengths 

The stability of a between base period lengths is 
measured by how much a varies in producing minimum figures 
of merit for equivalent validation years. The analysis will 
therefore observe optimal a for validation year 1 using a 
3-, 4-, and 5-year base period length, as well as validation 
year 2 using a 3- and 4-year base period length. 

Note that when reading the graphs for validation 
year one, the solid line is for the 3-year base, the the 
dot-dashed line is for the 4-year base and the dashed line 
is for the 5-year base. The graphs for validation year two 
have a solid line for the 3-year base and a dashed line for 
the 4-year base. The year corresponding to each type of 
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Figure 3. 13 Between VY One Stability -- Aviation 1LTS 
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Figure 3. 14 Between VY Two Stability -- Aviation 1LTS 
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Figure 3. 15 Between VY One Stability 



Aviation LTCOLS 
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Figure 3. 16 Between VY Two Stability — Aviation LTCOLS 
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Figure 3. 18 Between VY Two Stability 



Combat Spt. 1LTS 
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Figure 3. 17 Between VY One Stability -- Combat Spt. 1LTS 
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Figure 3. 19 Between VY One Stability -- Combat Spt. LTCOLS 
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Figure 3.20 Between VY Two Stability 



Combat Spt. LTCOLS 
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Figure 3.21 Between VY One Stability -- Ground Combat 1LTS 
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Figure 3.22 Between VY Two Stability -- Ground Combat 1LTS 
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Figure 3.23 Between VY One Stability -- Ground Combat LTCOLS 
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Figure 3.24 Between VY Two Stability -- Ground Combat LTCOLS 
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line is seen on one of the graphs, usually the 'MLE' graph, 
on each page. 

A review of Figures 3.13-3.24 show mixed results for 
the stability of the smoothing constants. The aviation 
aggregate (see Figures 3.13, 3.14, 3.15, and 3.16) shows a 
large instability between base period lengths for validation 
year one and an almost perfect stability for validation year 
two. This is again explained by the initiation of the 
Aviation Officer Continuation Pay program in 1981. The 
instability for validation year one is caused by the large 
decrease in optimal a for the 5-year base period from that 
of the other two. As discussed in the previous section, the 
5-year base contains 1981 data, so the « value no longer has 
to be 1. 0 to produce optimal figures of merit as it did in 
the 3- and 4-year bases (which explains the stability of a 
for validation year two). Thus, the instability of a is 
welcome here, as it is indicating a change in the aviation 
attrition environment. 

From Figures 3. 17 and 3. 18, we see that a is stable 
for combat support 1LTS. The only exception is the 
James-Stein estimator for validation year two. This excep- 
tion has a percent differential of 10%, while all of the 
other cases have percent differentials well below the arbi- 
trary acceptable level of 5%. Figures 3. 19 and 3. 20 show 
that the opposite is true for combat support LTCOLS; only 
the MLE shows stability. The percent differentials for TSCA, 
JS, and LTJS are all above 10% for the validation year one 
case, and have surprisingly small, though unacceptable, 
values of 5. 1%, 6. 7%, and 5. 6%, respectively, for the vali- 
dation year two case. They are "surprisingly" small because 
the ranges of optimal a are very large, .28, .36, and .32, 
respectively, and one would think that the percent differen- 
tial would be much bigger in light of these large ranges. A 
pairwise analysis of the optimal a values leads us to 
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believe that the optimal « for the 3-year base is signifi- 
cantly smaller than that for the 4- and 5-year bases, for 
which a appears stable. Thus, it would seem from this data- 
limited analysis that the optimal a values for TSCA, JS, and 
LTJS for combat support LTCOLS is stabilizing as the later 
years are included in the base. 

Figures 3. 21 and 3. 22 show a very stable optimal a 
value for all estimators in both validation year cases for 
ground combat 1LTS, except TSCA for validation year one, 
which has a percent differential of 8. 1%. As mentioned in 
the previous section, the reasons for the preponderance of a 
values of 1. 0 for this study block are unknown. All that is 
known is that these large a values indicate a very turbulent 
environment for the attrition of ground combat 1LTS, with 
the patterns changing dramatically from year to year. An 
analysis of Figures' 3. 23 and 3.24 shows that a is stable for 
ground combat LTCOLS as well. The only instability seen in 
these graphs is the MLE for validation year one, which has a 
percent differential of 6. 5%. Therefore, the conclusion is 
that a is stable for the ground combat aggregate. 
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IV. RESULTS 



A. GENERAL 

This chapter displays various data results, using a 
three-year base period, for the six estimators used in this 
study. 

B. COMPARISON WITH RESULTS FROM EARLIER WORKS 

Table 4 thru Table 6 display the figures of merit 
obtained by using a 3-year base period and smoothing the 
empirical rate of year 4 ( 1980) onto it. The figures of 
merit and optimal a represent those corresponding to the 
minimum transformed figures of merit for validation year 2 
(1982), since it was determined earlier in this study that 
the ability to forecast two years into the future should be 
the major concern of the Marine Corps. Table 7 thru Table 9 
show the results obtained by Robinson in his thesis [ Ref. 2: 
pp. 39-41] . These estimates are the same as those of a 
4-year base period with an « value of 0. In the analysis 
that follows, this estimation scheme is referred to as the 
"Robinson method" or the "Robinson estimation scheme. " 

One may notice a slight difference in the original scale 
figures of merit between Tables 7, 8, and 9 and Tables 5, 6, 
and 7 in Robinson's thesis [Ref. 2: pp. 39-41]. This is 
because a small error was found and corrected in APL func- 
tion RISKO since the submission of Robinson's thesis in 
March, 1986 ( that being the necessary addition of variable 
NV on lines 17-18 of the new RISKO, seen in Figure A. 8, 
which was absent from Robinson's version) [Ref. 2]. 
However, the errors resulting from this mistake in RISKO in 
Robinson's original scale figures of merit are quite small, 
and in no way invalidate the comparisons he made [Ref. 2]. 
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The results of using exponential smoothing are now 
presented for each of the six study blocks. The MOE for 
comparison is transformed FOM. 



EXPONENTIAL 


TABLE 4 

SMOOTHING AVIATION FIGURES 


OF MERIT 




a 


TRANSFORMED FOM 
1981 1982 


1983 


1st LT 


AGG ORIG 


i 


3. 976 


5. 634 


6. 260 


AGG TRANS 


0 


206. 877 


212. 181 


219. 619 


MLE 


1 


2. 437 


5. 746 


5. 156 


TSCA 


1 


1. 243 


5. 249 


4. 886 


JS 


1 


1. 402 


5. 360 


5. 126 


LTJS 


1 


1. 386 


5. 398 


5. 104 


LTCOL 


AGG ORIG 


1 


7. 898 


14. 829 


13. 238 


AGG TRANS 


1 


24. 733 


29. 123 


29. 083 


MLE 


1 


3. 423 


6. 942 


7. 604 


TSCA 


1 


3. 341 


7. 395 


8. 210 


JS 


1 


3. 302 


7. 382 


8. 167 


LTJS 


1 


3. 321 


7. 378 


8. 191 




a 


ORIGINAL FOM 
1981 1982 


1983 


1st LT 


AGG ORIG 


l 


40. 488 


12. 314 


50. 363 


AGG TRANS 


0 


384. 833 


306. 151 


886. 744 


MLE 


l 


13. 431 


18. 361 


42. 519 


TSCA 


l 


13. 428 


20. 271 


25. 000 


JS 


l 


16. 923 


22. 791 


26. 475 


LTJS 


l 


15. 797 


22. 423 


26. 032 


LTCOL 


AGG ORIG 


l 


74. 977 


113. 067 


50. 567 


AGG TRANS 


l 


63. 264 


92. 066 


76. 353 


MLE 


l 


20. 801 


65. 982 


32. 998 


TSCA 


l 


37. 928 


44. 331 


28. 742 


JS 


l 


29. 665 


42. 718 


27. 425 


LTJS 


l 


34. 297 


43. 152 


27. 835 



1. Aviation 

From a comparison of Tables 4 and 7, it is seen that 
exponential smoothing outperforms the Robinson estimation 
scheme in all cases and for all validation years except the 
transformed aggregate estimator. This exception is not 
important at all, since the figures of merit for this 
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TABLE 5 

EXPONENTIAL SMOOTHING COMBAT SUPPORT FIGURES OF MERIT 

TRANSFORMED FOM 





a 


1981 


1982 


1983 


1st LT 


AGG ORIG 


i 


1. 949 


1. 868 


1. 581 


AGG TRANS 


0 


4. 086 


4.580 


3. 995 


MLE 


. 56 


2. 402 


3. 093 


2. 420 


TSCA 


. 60 


1. 617 


2. 319 


1. 588 


JS 


. 90 


1. 688 


2. 257 


1. 559 


LTJS 


. 64 


1. 646 


2. 336 


1. 574 


LTCOL 


AGG ORIG 


0 


1. 532 


1. 889 


1. 956 


AGG TRANS 


1 


3. 372 


3. 649 


3. 757 


MLE 


. 38 


1. 860 


2. 520 


2. 754 


TSCA 


. 34 


1. 025 


1. 923 


2. 024 


JS 


. 26 


. 941 


1. 699 


1. 732 


LTJS 


. 30 


.971 1.782 

ORIGINAL FOM 


1. 861 




a 


1981 


1982 


1983 


1st LT 


AGG ORIG 


1 


88. 732 


62. 779 


50. 792 


AGG TRANS 


0 


234. 581 


136. 699 


96. 426 


MLE 


. 56 


1262. 038 


2810. 514 


3858. 857 


TSCA 


. 60 


75. 051 


86. 051 


65. 598 


JS 


. 90 


83. 608 


77. 961 


47. 631 


LTJS 


. 64 


75. 074 


85. 040 


60. 126 


LTCOL 


AGG ORIG 


0 


39. 172 


38. 241 


30. 186 


AGG TRANS 


1 


41. 704 


38. 615 


33. 363 


MLE 


. 38 


136. 467 


38. 763 


157. 663 


TSCA 


. 34 


23. 446 


31. 130 


29. 101 


JS 


. 26 


22. 863 


29. 503 


26. 596 


LTJS 


. 30 


23. 074 


29. 367 


27. 342 



estimator in all cases are so much higher than those for the 
others. This result extends into original scale, where 
exponential smoothing is again seen to have lower figures of 
merit. Special notice should be given to the exceptional 
performance of exponential smoothing in forecasting rates 
two and three years into the future as compared to the 
Robinson method. 

A major reason that exponential smoothing does so 
well in this case is the inflexibility of the Robinson esti- 
mation scheme. His scheme cannot readily adjust to changes 
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TABLE 6 

EXPONENTIAL SMOOTHING GROUND COMBAT FIGURES OF MERIT 

TRANSFORMED FOM 





a 


1981 


1982 


1983 


1st LT 


AGG ORIG 


. 92 


3. 045 


3. 663 


5. 988 


AGG TRANS 


0 


19. 842 


20.724 


24. Oil 


MLE 


. 94 


1. 858 


2. 165 


4. 555 


TSCA 


1 


1. 431 


1. 916 


4. 424 


JS 


1 


1. 472 


2. 008 


4. 566 


LTJS 


1 


1. 405 


1. 869 


4. 415 


LTCOL 


AGG ORIG 


1 


3. 540 


3. 427 


3. 671 


AGG TRANS 


0 


14. 393 


13. 453 


13. 777 


MLE 


. 42 


1. 389 


2. 293 


3. 085 


TSCA 


. 56 


1. 510 


2. 392 


3. 169 


JS 


. 58 


1. 505 


2. 313 


3. 039 


LTJS 


. 56 


1.494 2.358 

ORIGINAL FOM 


3. 108 




a 


1981 


1982 


1983 


1st LT 


AGG ORIG 


. 92 


64. 779 


62. 206 


80. 907 


AGG TRANS 


0 


235. 108 


236. 736 


303. 748 


MLE 


. 94 


175558. 747 


449. 816 


3491. 851 


TSCA 


1 


75. 514 


57. 155 


124. 972 


JS 


1 


82. 814 


55. 067 


122. 684 


LTJS 


1 


78. 815 


53. 724 


128. 000 


LTCOL 


AGG ORIG 


1 


90. 948 


110. 793 


69. 991 


AGG TRANS 


0 


221. 432 


197. 458 


128. 373 


MLE 


. 42 


37. 820 


1383. 609 


38. 934 


TSCA 


. 56 


40. 713 


52. 936 


49. 449 


JS 


. 58 


40. 615 


53. 345 


48. 090 


LTJS 


. 56 


40. Ill 


52. 986 


48. 151 



in the environment such as the Aviation Officer Continuation 
Pay program initiated in 1981. This program, as has been 
mentioned before, is believed to have radically changed 
attrition patterns for aviators as compared to 1977-1980. 
The exponential smoothing model's ability to anticipate this 
change, which it does by completely eliminating the effects 
of data from years 1977-1979 by having a values of 1, allows 
it to better predict attrition rates for the years following 
the change in the aviation environment. Eventually, once a 
base period of post-1981 data is established, the a values 
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TABLE 7 

FOUR-YEAR BASE ESTIMATE AVIATION FIGURES OF MERIT 

TRANSFORMED FOM 





1981 


1982 


1983 


1st LT 


AGG ORIG 


6. 405 


9. 645 


11. 393 


AGG TRANS 


197. 621 


202.022 


208. 808 


MLE 


3. 914 


9. 981 


10. 420 


TSCA 


3. 461 


9. 574 


10. 042 


JS 


3. 678 


9. 768 


10. 318 


LTJS 


3. 642 


9. 764 


10. 279 


LTCOL 


AGG ORIG 


9. 957 


17. 997 


15. 488 


AGG TRANS 


25. 093 


29. 506 


29. 310 


MLE 


4. 366 


9. 058 


8. 210 


TSCA 


5. 777 


10. 967 


10. 394 


JS 


5. 737 


10. 911 


10. 355 


LTJS 


5. 744 


10. 901 10. 340 

ORIGINAL FOM 




1981 


1982 


1983 


1st LT 


AGG ORIG 


31. 499 


34. 924 


54. 911 


AGG TRANS 


321. 894 


286. 068 


572. 830 


MLE 


22. 804 


46. 751 


57. 475 


TSCA 


19. 333 


38. 458 


49. 389 


JS 


21. 144 


39. 596 


51. 013 


LTJS 


20. 732 


39. 434 


50. 690 


LTCOL 


AGG ORIG 


74. 697 


111. 913 


51. 256 


AGG TRANS 


23050. 668 


248. 171 


3900. 284 


MLE 


33. 932 


55. 925 


22. 116 


TSCA 


37. 974 


57. 475 


37. 637 


JS 


36. 708 


54. 901 


35. 218 


LTJS 


37. 180 


55. 414 


35. 619 



should decrease substantially, thereby giving the base esti- 
mate some meaning in the estimation of the attrition rates. 
Further study into this matter should be undertaken as more 
data becomes available. 

2. Combat Support 

From a comparison of tables 5 and 8, one sees that 
in all cases for combat support 1LTS, the Robinson method 
figures of merit are slightly smaller than those produced by 
exponential smoothing. The values of a vary widely over the 
six estimators. Notice that maximum likelihood estimation 
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TABLE 8 



FOUR-YEAR BASE 


ESTIMATE COMBAT 
MERIT 


SUPPORT 


FIGURES OF 




TRANSFORMED 


FOM 




1981 


1982 


1983 


1st LT 


AGG ORIG 


1. 528 


1. 842 


1. 385 


AGG TRANS 


3. 308 


3. 722 


3. 238 


MLE 


2. 273 


2. 843 


2. 409 


TSCA 


1. 383 


2. 008 


1. 468 


JS 


1. 470 


2. 066 


1. 532 


LT JS 


1. 428 


2. 045 


1. 482 


LTCOL 


AGG ORIG 


1. 316 


1. 701 


1. 755 


AGG TRANS 


2. 601 


2. 898 


2. 981 


MLE 


1. 589 


2. 157 


2. 444 


TSCA 


. 910 


1. 689 


1. 752 


JS 


. 809 


1. 514 


1. 556 


LTJS 


. 831 


1. 560 


1. 611 




ORIGINAL FOM 




1981 


1982 


1983 


1st LT 


AGG ORIG 


79. 521 


69. 344 


54. 319 


AGG TRANS 


288. 130 


131. 044 


104. 613 


MLE 


141. 815 


145. 434 


104. 668 


TSCA 


73. 190 


85. 451 


71. 852 


JS 


74. 158 


81. 133 


63. 562 


LTJS 


74. 148 


84. 524 


68. 079 


LTCOL 


AGG ORIG 


34. 631 


46. 427 


38. 481 


AGG TRANS 


41. 785 


44. 739 


64. 345 


MLE 


36. 015 


62. 882 


48. 468 


TSCA 


27. 773 


40. 777 


45. 453 


JS 


26. 191 


36. 758 


41. 956 


LTJS 


26. 171 


36. 865 


42. 044 



produces alarmingly high figures of merit in the original 
scale. This is a direct result of the unsuitability of using 
the optimal a for transformed scale to predict original 
scale rates discussed in Chapter 3 (see also Appendix C). 
Additionally, its figures of merit in transformed scale are 
higher than those for all but one of the other estimators, 
so it appears that MLE is not a very good alternative to use 
for this study block. 
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TABLE 9 

FOUR-YEAR BASE ESTIMATE GROUND COMBAT FIGURES OF MERIT 



TRANSFORMED FOM 





1981 


1982 


1983 


1st LT 


AGG GRIG 


2. 774 


3. 609 


5. 720 


AGG TRANS 


18. 280 


19.064 


22. 044 


MLE 


2. 693 


4. 258 


6. 178 


TSCA 


1. 957 


3. 447 


5. 231 


JS 


2. 169 


3. 629 


5. 505 


LTJS 


2. 045 


3. 482 


5. 334 


LTCOL 


AGG ORIG 


3. 692 


3. 596 


3. 783 


AGG TRANS 


13. 108 


12. 209 


12. 546 


MLE 


1. 252 


2. 029 


2. 925 


TSCA 


1. 598 


2. 414 


3. 333 


JS 


1. 584 


2. 332 


3. 223 


LTJS 


1. 567 


2.367 3.263 

ORIGINAL FOM 




1981 


1982 


1983 


1st LT 


AGG ORIG 


79. 406 


87. 445 


110. 324 


AGG TRANS 


271. 834 


255. 215 


341. 258 


MLE 


91. 399 


134. 172 


170. 424 


TSCA 


69. 166 


94. 099 


118. 799 


JS 


76. 811 


98. 672 


125. 292 


LTJS 


73. 141 


94. 657 


121. 432 


LTCOL 


AGG ORIG 


100. 321 


122. 612 


75. 718 


AGG TRANS 


250. 966 


342. 375 


4494. 138 


MLE 


35. 833 


45. 817 


56. 710 


TSCA 


37. 521 


51. 760 


57. 158 


JS 


36. 909 


50. 712 


55. 615 


LTJS 


36. 373 


50. 828 


55. 562 



Comparing Tables 5 and 8 for lieutenant colonels, we 
see that the Robinson method's estimates are barely better 
than the exponential smoothing estimates in all cases. The 
values of a producing the minimum figures of merit for 1982 
range from .26 to . 38 over the 4 estimators introduced by 
Robinson [Ref. 2]. However, as one can see from Figure 3.4, 
the curves are flat enough in this range that choosing any a 
in this range would produce good results. 

Exponential smoothing, therefore, proves to be an 
good technique for predicting attrition rates for the combat 
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support aggregate. Its forecasting ability is almost as good 
as Robinson's method while requiring far less data. 

3. Ground Combat 

A comparison of Tables 6 and 9 shows that exponen- 
tial smoothing outperforms the Robinson method's estimates 
in most cases for the ground combat aggregate. For 1LTS, 
exponential smoothing is better in all cases except for the 
aggregate estimators. The improvement over the Robinson 
method figures of merit is almost 50% for validation year 
1982 for MLE, TSCA, JS, and LTJS, and is also quite notice- 
able for the other two validation years. Again, the reasons 
for the a values near or at 1. 0 are unknown, but they indi- 
cate a radical change in the attrition environment for 
ground combat 1LTS. Further study is needed to determine 
whether or not this turbulence is specific to our data. But 
whatever causes this apparent yearly change in attrition 
patterns, the exponential smoothing model shows its efficacy 
as a forecasting model in this study block by anticipating 
this change and discounting the now-irrelevant base period 
data, with the result being significant improvements in 
forecasting ability over the methods of Robinson. 

For ground combat LTCOLS, The AGG ORIG, TSCA, JS, 
and LTJS estimators using exponential smoothing have smaller 
figures of merit than the Robinson estimation scheme, and 
the figures of merit for AGG TRANS and MLE are slightly 
higher for exponential smoothing than for the Robinson 
method. However, the differences are small in both direc- 
tions; they are not nearly as pronounced for this study 
block as they are for combat support 1LTS. These results 
hold for all three validation years. Again, we see the expo- 
nential smoothing model producing attrition rates as good 
as, if not better than, a method which theoretically 
requires much more data. This study block, then, also shows 
exponential smoothing to be a good model to perform the 
attrition rate calculations. 
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C. ATTRITION RATES 



The attrition rates for all 7 operational fields for 
first lieutenants and lieutenant colonels are presented in 
Tables 10 thru 23. Attrition rates are computed for all 6 
estimators and all 31 lengths of service. 

Unlike the estimates produced by Robinson for the orig- 
inal aggregate [Ref. 2: pp. 54-67], those calculated by 
exponential smoothing are not the same for all cells with 
non-zero rates, except when «=0, as is the case with combat 
support LTCOLS. They would be the same were the smoothing 
done in original scale, but since it is done in transformed 
scale, variability is introduced. This is because taking the 
average of two sets of data in transformed scale, linearly 
combining them using Equation 2.1, and then transforming 
back to original scale will not produce a single average for 
all cells as would linearly combining two original scale 
average rates. This is basically the same reason that there 
is variability in the transformed aggregate estimates. 
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TABLE 10 



AVIATION ATTRITION RATES FOR 1ST LTS 



lGG orig 


AGG TRANS 


MLE 


TSCA 


JS 


LTJS 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


. 0179 


. 1154 


. 0008 


. 0003 


. 0005 


. 0005 


. 0253 


. 2283 


. 0266 


. 0283 


. 0300 


. 0294 


. 0159 


. 2527 


. 0090 


. 0100 


. 0114 


. 0106 


. 0282 


. 2221 


. 0402 


. 0421 


. 0437 


. 0434 


. 0216 


. 1442 


. 0240 


. 0270 


. 0278 


. 0278 


. 0204 


. 0758 


. 0386 


. 0432 


. 0430 


. 0430 


. 0350 


. 0064 


. 0621 


. 0729 


. 0698 


. 0698 


. 0788 


. 0067 


. 0251 


. 0533 


. 0510 


. 0510 


. 0710 


. 0072 


. 1197 


. 1335 


. 1282 


. 1282 


. 0480 


. 0182 


. 0063 


. 0307 


. 0277 


. 0277 


. 0123 


0 


. 0033 


. 0155 


. 0097 


. 0109 


. 0832 


0 


. 0284 


. 2679 


. 2027 


. 2313 


. 0725 


0 


. 5000 


. 5000 


. 5000 


. 5000 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 
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TABLE 11 

AVIATION ATTRITION RATES FOR LTCOLS 



AGG ORIG 


AGG TRANS 


MLE 


TSCA 


JS 


LTJS 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


. 0098 


0 


. 2957 


. 0298 


. 0373 


. 0373 


. 0089 


. 1001 


. 0001 


. 0101 


. 0076 


. 0076 


. 0907 


. 0672 


. 0270 


. 0010 


. 0003 


. 0004 


. 1304 


. 1762 


. 0015 


. 0026 


. 0046 


. 0034 


. 1302 


. 2278 


. 0015 


. 0012 


. 0029 


. 0016 


. 1442 


. 2514 


. 0717 


. 0790 


. 0850 


. 0820 


. 1739 


. 2575 


. 1568 


. 1642 


. 1684 


. 1681 


. 1791 


. 2639 


. 1991 


. 2054 


. 2084 


. 2084 


. 1812 


. 2512 


. 2447 


. 2506 


. 2514 


. 2514 


. 1919 


. 2182 


. 2464 


. 2549 


. 2542 


. 2542 


. 1776 


. 1835 


. 2676 


. 2765 


. 2735 


. 2735 


. 1281 


. 1071 


. 3326 


. 3401 


. 3303 


. 3313 


. 1018 


. 0237 


. 2925 


. 3070 


. 2926 


. 2950 


. 1602 


. 0016 


. 6575 


. 6311 


. 5953 


. 6145 


. 1789 


. 0430 


. 4119 


. 4343 


. 3998 


. 4139 


. 1184 


. 0500 


. 9944 


. 9683 


. 9423 


. 9605 



75 



LOS 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 



TABLE 12 



COMBAT SUPPORT ATTRITION RATES FOR 1ST LTS 
CODE 07 ENGINEERS 



AGG ORIG 


AGG TRANS 


MLE 


TSCA 


JS 


LTJS 


0 


0 


0 


0- 


0 


0 


0 


0 


0 


0 


0 


0 


. 2357 


. 3571 


. 0901 


. 1119 


. 1388 


. 1134 


. 1987 


. 3973 


. 2610 


. 2697 


. 3124 


. 2791 


. 1547 


. 3940 


. 1446 


. 1455 


. 1768 


. 1454 


. 1903 


. 3793 


. 1838 


. 1861 


. 2144 


. 1897 


. 1870 


. 3511 


. 2990 


. 3072 


. 3086 


. 3135 


. 2066 


. 2976 


. 3226 


. 3443 


. 3214 


. 3325 


. 1820 


. 1756 


. 3753 


. 4145 


. 2619 


. 3768 


. 1824 


. 1916 


. 0520 


. 1312 


. 0871 


. 1266 


. 1933 


. 2023 


. 2075 


. 2523 


. 2371 


. 2444 


. 2578 


. 1953 


. 2667 


. 3089 


. 3358 


. 3200 


. 2370 


. 1711 


. 0411 


. 1327 


. 1338 


. 1326 


. 1915 


. 1461 


0 


. 0801 


. 0843 


. 0868 


. 1216 


. 1711 


. 0064 


. 0269 


. 0273 


. 0329 


. 1375 


. 1618 


. 0036 


. 0373 


. 0384 


. 0463 


. 1683 


. 1218 


. 0005 


. 0748 


. 0640 


. 0754 


. 1616 


. 1149 


. 0009 


. 0739 


. 0586 


. 0723 


. 2148 


. 1461 


. 0012 


. 0966 


. 1065 


. 0977 


. 2209 


. 1218 


. 0018 


. 1161 


. 1132 


. 1113 


. 2370 


. 1002 


. 0646 


. 1989 


. 1384 


. 1692 


. 2869 


. 1077 


. 0924 


. 2503 


. 1946 


. 2252 


. 2756 


. 0923 


. 0118 


. 1912 


. 1728 


. 1845 


. 2062 


. 0661 


. 0006 


. 1463 


. 1005 


. 1220 


. 2815 


. 0168 


. 1583 


. 3626 


. 4550 


. 3571 


. 2370 


. 0003 


. 1671 


. 4058 


. 4575 


. 3916 


. 2370 


. 0003 


. 1671 


. 4058 


. 4575 


. 3916 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 
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TABLE 13 



COMBAT SUPPORT ATTRITION RATES FOR 1ST LTS 
CODE 13 COMMUNICATIONS 



AGG ORIG 


AGG TRANS 


MLE 


TSCA 


JS 


LTJS 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


. 2370 


. 3579 


. 1800 


. 1738 


. 2067 


. 1802 


. 1939 


. 4065 


. 2699 


. 2669 


. 2949 


. 2719 


. 1855 


. 4031 


. 2084 


. 2063 


. 2264 


. 2063 


. 2148 


. 3893 


. 3155 


. 3184 


. 3204 


. 3245 


. 2086 


. 3540 


. 2932 


. 3072 


. 3356 


. 3087 


. 2077 


. 2932 


. 4424 


. 4523 


. 3931 


. 4370 


. 2561 


. 2847 


. 2528 


. 2799 


. 2783 


. 2767 


. 2291 


. 3230 


. 1881 


. 2078 


. 2178 


. 2172 


. 2091 


. 3304 


. 1381 


. 1514 


. 2213 


. 1667 


. 1977 


. 2909 


. 0941 


. 1268 


. 1807 


. 1422 


. 1715 


. 2179 


. 0097 


. 0579 


. 0681 


. 0663 


. 1332 


. 1077 


. 0175 


. 0942 


. 0410 


. 0857 


. 1350 


. 0367 


. 0040 


. 1113 


. 0398 


. 0774 


. 0905 


. 0003 


. 0265 


. 1909 


. 0222 


. 1426 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


. 2370 


. 0003 


. 1671 


. 4058 


. 4575 


. 3916 


. 2370 


. 0003 


. 1671 


. 4058 


. 4575 


. 3916 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 
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0 


0 


0 
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TABLE 14 



COMBAT SUPPORT ATTRITION RATES FOR 1ST LTS 
CODE 20 MOTOR TRANSPORT 



GG ORIG 


AGG TRANS 


MLE 


TSCA 


JS 


LTJS 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


. 0905 


. 1618 


. 0140 


. 0173 


. 0102 


. 0172 


. 0983 


. 2472 


. 2572 


. 2942 


. 2035 


. 2776 


. 0926 


. 2690 


. 2669 


. 2635 


. 2461 


. 2534 


. 0812 


. 2751 


. 1576 


. 2076 


. 2071 


. 2178 


. 0670 


.2149 


. 2195 


. 2537 


. 1756 


. 2345 


. 0817 


. 0752 


. 1060 


. 2807 


. 1599 


. 2411 


. 1616 


. 0266 


. 9991 


. 8392 


. 5214 


. 8041 


. 1271 


. 0566 


. 9236 


. 8146 


. 6270 


. 7987 


. 1746 


. 0367 


. 0002 


. 1614 


. 0738 


. 1243 


. 1721 


. 1002 


. 1863 


. 3277 


. 2873 


. 3184 


. 1721 


. 1002 


. 0003 


. 0938 


. 0684 


. 0805 


. 1528 


. 1461 


. 0017 


. 0544 


. 0508 


. 0652 


. 1862 


. 1404 


0 


. 0839 


. 0803 


. 0881 


. 2692 


. 1798 


. 0102 


. 1228 


. 1626 


. 1249 


. 2043 


. 1345 


. 0005 


. 0949 


. 0966 


. 0950 


. 1862 


. 0840 


0 


. 1137 


. 0811 


. 0938 


. 2546 


. 0661 


. 0070 


. 1956 


. 1508 


. 1759 


. 2148 


. 0752 


. 0012 


. 1508 


. 1094 


. 1288 


. 3091 


. 0467 


. 3740 


. 4019 


. 4646 


. 3950 


. 2546 


. 0661 


. 0070 


. 1956 


. 1508 


. 1759 


. 3011 


. 0367 


. 1552 


. 3486 


. 4546 


. 3463 


. 3011 


. 0367 


. 1552 


. 3233 


. 4502 


. 3243 


. 2692 


. 0081 


. 1605 


. 3767 


. 4561 


. 3686 


. 2692 


. 0081 


. 1605 


. 3731 


. 4555 


. 3654 


. 2692 


. 0081 


. 1605 


. 3731 


. 4555 


. 3654 


. 1683 


. 0018 


. 0005 


. 2301 


. 0703 


. 1769 


. 1216 


. 0088 


. 0516 


. 2748 


. 0477 


. 2210 


. 2148 


. 0088 


. 1726 


. 4482 


. 4620 


. 4272 


. 2692 


. 0081 


. 1605 


. 3559 


. 4525 


. 3500 
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TABLE 15 



COMBAT SUPPORT ATTRITION RATES FOR LTCOLS 





CODE 07 


ENGINEERS 






AGG ORIG 


AGG TRANS 


MLE 


TSCA 


JS 


LTJS 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


. 2019 


. 0094 


. 1263 


. 3156 


. 3495 


. 3273 


0 


0 


0 


0 


0 


0 


. 1404 


. 1507 


. 0001 


. 1012 


. 1084 


. 1100 


. 1531 


• .2705 


. 0001 


. 0370 


. 0780 


. 0546 


. 1542 


. 2939 


. 0001 


. 0299 


. 0739 


. 0443 


. 1422 


. 3012 


. 2587 


. 2506 


. 2557 


. 2589 


. 1328 


. 2984 


. 0443 


. 0897 


. 1303 


. 1100 


. 1408 


. 2908 


. 1488 


. 1486 


. 1722 


. 1660 


. 1514 


. 2662 


. 0962 


. 1260 


. 1492 


. 1452 


. 1366 


. 2368 


. 2638 


. 2742 


. 2543 


. 2591 


. 1492 


. 1979 


. 1442 


. 2581 


. 2640 


. 2534 


. 1746 


. 0977 


. 0007 


. 1273 


. 1289 


. 1237 


. 1475 


. 0866 


. 6841 


. 6243 


. 5213 


. 5931 


. 1289 


. 0218 


. 6117 


. 5407 


. 4169 


. 4692 


. 1289 


. 0218 


. 0685 


. 3187 


. 2386 


. 2514 


. 1206 


. 0748 


. 5488 


. 4937 


. 4012 


. 4345 
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TABLE 16 



COMBAT SUPPORT ATTRITION RATES FOR LTCOLS 
CODE 13 COMMUNICATIONS 



AGG ORIG 


AGG TRANS 


MLE 


TSCA 


JS 


LT JS 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


Q 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


. 1385 


. 0010 


. 0777 


. 4148 


. 2942 


. 3253 


. 1200 


. 1767 


. 0052 


. 1298 


. 1252 


. 1324 


. 1304 


. 2615 


. 0012 


. 0530 


. 0831 


. 0681 


. 1535 


. 3052 


. 0001 


. 0273 


. 0723 


. 0403 


. 1381 


. 3222 


. 1480 


. 1784 


. 2087 


. 2010 


. 1238 


. 3170 


. 1689 


. 1914 


. 2051 


. 1967 


. 1515 


. 3170 


. 1590 


. 1927 


. 2193 


. 2170 


. 1673 


. 2841 


. 1437 


. 1654 


. 1799 


. 1763 


. 1495 


. 2615 


. 2642 


. 2531 


. 2573 


. 2559 


. 1488 


. 2302 


. 1832 


. 2521 


. 2699 


. 2524 


. 1369 


. 1644 


. 3054 


. 3026 


. 2411 


. 2668 


. 1637 


. 1177 


. 3865 


. 4655 


. 3164 


. 3906 


. 1475 


. 0866 


0 


. 1736 


. 1510 


. 1499 


. 1475 


. 0866 


. 1387 


. 3213 


. 2692 


. 2749 


. 1321 


. 2368 


. 3756 


. 3945 


. 3490 


. 3531 
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TABLE 17 



COMBAT SUPPORT ATTRITION RATES FOR LTCOLS 
CODE 20 MOTOR TRANSPORT 



AGG ORIG 


AGG TRANS 


MLE 


TSCA 


JS 


LTJS 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


. 1780 


. 0010 


. 1573 


. 3692 


. 3769 


. 3584 


. 1780 


. 0010 


. 1573 


. 3692 


. 3769 


. 3584 


. 1637 


. 0094 


. 0003 


. 3062 


. 2117 


. 2227 


. 1475 


. 0354 


0 


. 2113 


. 1550 


. 1551 


. 1385 


. 0623 


. 1923 


. 3483 


. 2709 


. 2801 


. 1238 


. 0491 


. 4323 


. 4363 


. 3445 


. 3686 


. 1424 


. 0491 


. 2332 


. 4020 


. 3082 


. 3311 


. 1780 


. 0010 


. 1573 


. 3692 


. 3769 


. 3584 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


. 1637 


. 0094 


. 0003 


. 3062 


.2117 


. 2227 


. 1570 


. 0623 


. 8041 


. 6442 


. 4519 


. 5405 
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TABLE 18 



GROUND COMBAT ATTRITION RATES FOR 1ST LTS 
CODE 03 INFANTRY 



AGG ORIG AGG TRANS 
0 0 



. 1835 


0 


. 1861 


. 3367 


. 1690 


. 3828 


. 1449 


. 3829 


. 1808 


. 3725 


. 1792 


. 3453 


. 1474 


. 2996 


. 1859 


. 2312 


. 2204 


. 2653 


. 1974 


. 2649 


. 1847 


. 2271 


. 1927 


. 1047 


. 1722 


. 0208 


. 1913 


. 0566 


. 2367 


. 0566 


. 1949 


. 1904 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 



MLE 


TSCA 


JS 


LT JS 


0 


0- 


0 


0 


. 4514 


. 5000 


. 5000 


. 5000 


. 0129 


. 0217 


. 0385 


. 0246 


. 1093 


. 1063 


. 1291 


. 1106 


. 1458 


. 1423 


. 1632 


. 1472 


. 1732 


. 1764 


. 1937 


. 1822 


. 2387 


. 2407 


. 2486 


. 2486 


. 2559 


. 2637 


. 2630 


. 2630 


. 0855 


. 1026 


. 1105 


. 1105 


. 2155 


. 2353 


. 2333 


. 2333 


. 1681 


. 1832 


. 1870 


. 1870 


. 0814 


. 1031 


. 1105 


. 1105 


. 1769 


. 2217 


. 1971 


. 2000 


. 0038 


. 0694 


. 0547 


. 0547 


. 0094 


. 1855 


. 1114 


. 1414 


. 0570 


. 2892 


. 1882 


. 2366 


. 0113 


. 2494 


. 1371 


. 1898 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 
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TABLE 19 



GROUND COMBAT ATTRITION RATES FOR 1ST LTS 
CODE 05 ARTILLERY 



iGG ORIG 


AGG TRANS 


MLE 


TSCA 


JS 


LTJS 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


. 1879 


. 2913 


. 0126 


. 0285 


. 0435 


. 0328 


. 1578 


. 3427 


. 1780 


. 1797 


. 1932 


. 1869 


. 1520 


. 3442 


. 2301 


. 2318 


. 2404 


. 2396 


. 1682 


. 3218 


. 1889 


. 1938 


. 2035 


. 2022 


. 1597 


. 2750 


. 2756 


. 2793 


. 2734 


. 2734 


. 1467 


. 1705 


. 4580 


. 4711 


. 4242 


. 4502 


. 1707 


. 1329 


. 1988 


. 2302 


. 2093 


. 2103 


. 1684 


. 1424 


. 2686 


. 2891 


. 2614 


. 2682 


. 1889 


. 1537 


. 1291 


. 1644 


. 1559 


. 1559 


. 1741 


. 0922 


. 1440 


. 1918 


. 1695 


. 1703 


. 1364 


. 0369 


. 1185 


. 1753 


. 1441 


. 1498 


. 0906 


. 0027 


. 1420 


. 2080 


. 1458 


. 1698 


. 1328 


. 2399 


. 0047 


. 1210 


. 0486 


. 0761 


. 2050 


0 


. 4504 


. 5000 


. 5000 


. 5000 


. 2050 


0 


. 4504 


. 5000 


. 5000 


. 5000 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 
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TABLE 20 



GROUND COMBAT ATTRITION RATES FOR 1ST LTS 
CODE 10 TANKS 



AGG ORIG 

0 

0 

. 2003 
. 1756 
. 1287 
. 1807 
. 1663 
. 1307 
. 1471 
. 2572 
. 2088 
. 1391 
. 1479 
. 1192 
. 0958 
0 
0 

. 2050 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 



AGG TRANS 

0 

0 

. 1607 
. 2534 
. 2 543 
. 2516 
. 2019 
. 1065 
. 0208 
. 0392 
. 0683 
. 0230 
. 0731 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 



MLE 


TSCA 


0 


o- 


0 


0 


. 0256 


. 0654 


. 1063 


. 1121 


. 1321 


. 1326 


. 2540 


. 2716 


. 2822 


. 2899 


. 0716 


. 0850 


0 


. 0355 


. 0830 


.2134 


. 0271 


. 1059 


. 0011 


. 0258 


. 0006 


. 1037 


. 0124 


. 1062 


. 0219 


. 1489 


0 


0 


0 


0 


. 4504 


. 5000 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 



JS 


LTJS 


0 


0 


0 


0 


0694 


. 0694 


1215 


. 1215 


1403 


. 1403 


2633 


. 2633 


2719 


. 2731 


0808 


. 0808 


0279 


. 0279 


1760 


. 1861 


0932 


. 0932 


0206 


. 0206 


0526 


. 0691 


0364 


. 0622 


0331 


. 0813 


0 


0 


0 


0 


5000 


. 5000 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 
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TABLE 22 



GROUND COMBAT ATTRITION RATES FOR LTCOLS 
CODE 05 ARTILLERY 



AGG ORIG 


AGG TRANS 


MLE 


TSCA 


JS 


LTJS 


0 


0 


0 


o- 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


. 0293 


. 0310 


. 0033 


. 0655 


. 0508 


. 0587 


. 1459 


. 1167 


0 


. 0336 


. 0362 


. 0361 


. 1415 


. 2141 


0 


. 0152 


. 0192 


. 0171 


. 1232 


. 2424 


. 0136 


. 0241 


. 0276 


. 0262 


. 1475 


. 2390 


. 2642 


. 2813 


. 2812 


. 2794 


. 1457 


. 2282 


. 0751 


. 0990 


. 1051 


. 1024 


. 1276 


. 2123 


. 1129 


. 1229 


. 1282 


. 1257 


. 1361 


. 1742 


. 1583 


. 1996 


. 2008 


. 1984 


. 1497 


. 1446 


. 3989 


. 4354 


. 4264 


. 4267 


. 1374 


. 0897 


. 2138 


. 2547 


. 2495 


. 2472 


. 1187 


. 0519 


. 4907 


. 4585 


. 4331 


. 4451 


. 0776 


. 0448 


. 8269 


. 7777 


. 7420 


. 7573 


. 0877 


. 0310 


. 2152 


. 2266 


. 1932 


. 2102 


. 2426 


. 1625 


. 6976 


. 5877 


. 5592 


. 5751 


. 2966 


. 0207 


. 6503 


. 5945 


. 5703 


. 5848 
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TABLE 21 



GROUND COMBAT ATTRITION RATES FOR LTCOLS 
CODE 03 INFANTRY 



AGG ORIG 


AGG TRANS 


MLE 


TSCA 


JS 


LTJS 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


. 0151 


. 1625 


. 0596 


. 1020 


. 0867 


. 0964 


. 0277 


. 0191 


. 0034 


. 0150 


. 0122 


. 0151 


. 1500 


. 2018 


. 0001 


. 0204 


. 0246 


. 0226 


. 1327 


. 3044 


0 


. 0054 


. 0089 


. 0061 


. 1195 


. 3290 


. 0029 


. 0058 


. 0092 


. 0064 


. 1311 


. 3383 


. 1187 


. 1211 


. 1281 


. 1237 


. 1476 


. 3345 


. 0848 


. 1002 


. 1091 


. 1027 


. 1411 


. 3277 


. 1395 


. 1592 


. 1673 


. 1623 


. 1314 


. 3099 


. 1477 


. 1667 


. 1736 


. 1702 


. 1627 


. 2941 


. 2292 


. 2490 


. 2528 


. 2508 


. 1466 


. 2590 


. 2207 


. 2388 


. 2410 


. 2396 


. 1394 


. 2159 


. 2107 


. 2034 


. 2009 


. 2039 


. 1345 


. 1618 


. 2265 


. 2134 


. 2062 


. 2114 


. 1707 


. 0897 


. 2299 


. 2495 


. 2424 


. 2415 


. 1741 


. 0661 


. 1263 


. 1549 


. 1471 


. 1504 


. 1666 


. 1463 


. 4337 


. 4319 


. 4168 


. 4225 
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TABLE 23 



GROUND COMBAT ATTRITION RATES FOR LTCOLS 
CODE 10 TANKS 



AGG ORIG 


AGG TRANS 


MLE 


TSCA 


JS 


LTJS 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


. 0972 


0 


. 0006 


. 2593 


. 2115 


. 2288 


. 1147 


. 0077 


. 0002 


. 1060 


. 0935 


. 0951 


. 1538 


. 0085 


. 0001 


. 0790 


. 0752 


. 0746 


. 1315 


. 0634 


0 


. 0396 


. 0402 


. 0404 


. 1032 


. 0661 


. 1689 


. 1690 


. 1629 


. 1637 


. 1956 


. 0548 


. 0596 


. 1936 


. 1939 


. 1888 


. 1769 


. 0689 


. 2050 


. 2508 


. 2441 


. 2422 


. 1716 


. 0489 


. 0317 


. 0899 


. 0879 


. 0878 


. 1205 


. 0429 


. 1361 


. 1596 


. 1460 


. 1556 


. 1323 


. 0085 


. 3439 


. 4017 


. 3756 


. 3854 


. 1995 


. 0017 


. 0889 


. 2828 


. 2869 


. 2775 


. 1850 


. 0077 


. 1695 


. 2636 


. 2373 


. 2452 


. 1083 


. 0448 


. 1714 


. 4064 


. 3781 


. 3828 


. 1252 


. 0876 


. 3967 


. 5706 


. 5365 


. 5446 


. 2061 


. 0004 


. 6177 


. 5570 


. 5312 


. 5376 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

This study investigated the ability of the exponential 
smoothing model to update attrition rates for Marine Corps 
manpower models from year to year. Its value as a yearly 
updating scheme has been demonstrated in this work. 

For the reader who wants a single general purpose value 
of a, we offer the value a = . 4, but do so reluctantly. 
Brown has suggested using an a between . 01 and . 3 in the 
applications he studied [ Ref. 6: p. 106] . Our work, 
involving manpower attrition rates, appears to call for 
larger values. However, we have identified some special 
situations for which the smoothing constant should be 
considerably larger, e. g. , aviation LTCOLS and 1LTS with a 
pre-1981 base forecasting post-1981 rates. 

Exponential smoothing produced estimates for the combat 
support and ground combat aggregates that were, more often 
than not, better than those produced by the methods of Major 
Robinson [Ref. 2], without needing, in theory, the massive 
data files his methods use to produce the attrition rate 
estimates. The exponential smoothing model reflected the 
change in the aviation environment that occurred in 1981 and 
easily outperformed the Robinson method's estimates in this 
aggregate for the years following because of the inflexi- 
bility of the Robinson method. It also anticipated an 
unknown source of turbulence in the ground combat first 
lieutenant attrition rates, and bested the estimates of 
Robinson's methods for this study block. 

Also seen in this study was that three of the four esti- 
mators presented by Robinson [Ref. 2], transformed cell 
scale average, James-Stein, and limited translation 
James-Stein, outperform the current method used by the 
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Marine Corps when exponential smoothing is used in all vali- 
dation years and study blocks except validation years 2 and 
3 for combat support 1LTS. These two exceptions, however, 
were also the only ones seen in Robinson's thesis [Ref. 2: 
p. 40] . Unfortunately, none of these three emerges as the 
clear-cut "best" estimator in this study, but all are better 
than the current Marine Corps estimator.- 

The maximum likelihood estimator shows the same good 
performance in transformed scale in this study as it did in 
Robinson's [Ref. 2: pp. 39-41]. However, there are more 
cases of large original scale figures of merit for MLE in 
the present work than were seen previously (see Tables 4-9). 
These extraordinarily large figures are seen for the ground 
combat aggregate for both grades and for combat support 
1LTS. Not coincidentally, these are the same study blocks 
for which the use of the optimal a for transformed scale to 
produce original scale estimates was determined to be unsui- 
table in the analysis of Appendix C. The ability of the MLE 
to produce attrition rates better than those of aggregate 
estimators for the other three study blocks indicates that 
it, too, may be better than the aggregate methods currently 
in use, and warrants further study into the cause of the 
problems in original scale just mentioned. Perhaps if 
smoothing is done in original scale rather than transformed 
scale, the optimal properties of MLE will be better 
displayed. 

Thus, Robinson's conclusion [Ref. 2: p. 68] that the 
MLE, TSCA, JS, and LTJS estimators should be given serious 
consideration for replacing the Marine Corps' current scheme 
is reiterated in this work. 

B. RECOMMENDATIONS 

Despite the encouraging results produced herein, it is 
not recommended at this time that either exponential 
smoothing nor any of the four promising estimators presented 
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as alternatives to the current method be implemented. 
Further study is needed in the following areas: 

1. Base Period. The problem of anticipating major changes 
in Marine Corps policy such as the one that the avia- 
tion OF underwent in 1981 and their effects on the 
length of the base period should be studied. Also, as 
more data becomes available, more analyses like the 
ones conducted herein to determine optimal base period 
length should be conducted to produce stronger conclu- 
sions about base period length. 

2. Aviation. As more data is obtained, the aviation 
aggregate should be re-examined and 1981 used as the 
first base year to see if lower values of a can be 
produced, indicating a consistency from year to year 
in the loss rates which did not exist from 1977 to 
1983. 

3. Ground Combat 1LTS. The reasons behind the values of a 
being 1.0, which indicate a year-to-year change in 
attrition patterns for this study block, should be 
investigated. Despite the excellent results achieved 
with respect to the low figures of merit produced, a 

E ersistence of such a values would indicate that the 
ase period is totally unimportant, thereby making the 
use of exponential smoothing unnecessary. The issue 
is, therefore, whether or not the tt values seen for 
ground combat 1LTS herein are the product of a 
peculiar set of data. A new analysis of this study 
block should be undertaken using years other than 
1977-1979 as the base period as soon as enough data 
becomes available to perform validations on newly 
produced attrition rates. 

4. Maximum Likelihood Estimation. In light of the prom- 
ising performance seen for this estimator in trans- 
formed scale using exponential smoothing, an 
investigation of the unsuitability of using trans- 
formed optimal « values to produce original scale loss 
rates seen for certain study blocks is needed. 

5. Aggregation. As recommended by Major Robinson [Ref. 2: 
p. 69j, the work of Amin Elseramegy [Ref. 5] , who used 
the CART routine to find aggregations with encouraging 
results, should be investigated. Also, the work of 
Major Tucker [Ref. 1: pp. 75-84], which demonstrated 

increased attrition at certain lengths of service for 
certain grades, should be followed up by attempting to 
aggregate by LOS instead of OF to see if the results 
can be improved. 
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APPENDIX A 

FUNCTIONS USED IN CALCULATIONS 



The following APL functions, and ones similar to them, 
are used in this study to produce the figures of merit and 
attrition rates seen in the tables in Chapters III and IV. 
These procedures require that the following variables be 
global, i. e. , defined throughout the workspace: 

(1) N = the estimation period central inventory array, 
e. g. , NA5, 

(2) Y = the estimation period loss array, e. g. , YA5, 

(3) VN = the validation period centarl inventory array, 
e. g. , VNA5, 

(4) VY = the validation period loss array, e. g. , VYA5, 

(5) AN = the estimation period average inventory array, 
e. g. , ANA5 , 

(6) G = the forced James-Stein shrinkage rate, and 

(7) DEE = the factor used in limited translation 
James-Stein shrinkage. 

Functions like A5B seen herein set the values of all of 
these parameters except G, which is set in the workspace 
itself to 0. 

In order to create the arrays of figures of merit for a 
values between 0 and 1 for a three-year base period, APL 
function ALPHAHAT is called (see Figure A. 1). ALPHAHAT calls 
A5B in Figure A. 2 to produce the 3-year base estimates from 
years 1977-1979 and the empirical estimates for 1980. To do 
this, A5B sets the values of the global variables, then 
calls ESTIM, seen in Figure A. 3, to produce the base esti- 
mates. ESTIM, in turn, calls BINPREP in Figure A. 4, SUMSQ in 
Figure A. 5, and MLE in Figure A. 6 in performing its calcula- 
tions. A5B then calls XFOUR in Figure A. 7 which resets the 
global variables N, Y, and AN, and calls ESTIM to produce 
the empirical transformed attrition figures. Notice that the 
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V ALPHAHAT : LB :UB: STEP; RATES ;M 1 ;M2 ; ALPHA 1 NUMSTEPS 
n THIS FUNCTION WILL FIND THE VALUES OF THE RISK 
r FUNCTION IN BOTH TRANSFORMED AND ORIGINAL SPACE 
r FOR VARYING VALUES OF ALPHA BETWEEN ANY TWO 
r LIMITS SET BY THE USER. 

R ' INPUT LOWER BOUND ' 

LB+ 0 

r' INPUT UPPER BOUND ' 

UB+ 1 

o ' INPUT STEPSIZE ' 

STEP* 0.02 
ALPHA+LB 

NUMSTEPS+1+ (UB-LB )*STEP 
r INITIALIZE ARRAYS ORl-ORL AND TR1-TR6, WHICH 
R ARE THE ARRAYS OF FOM'S FOR ALPHAS BETWEEN 
fl LB AND UB BY STEPSIZE STEP. VALUES 1-6 
r CORRESPOND TO AGG ORIG, AGG TRANS, MLE , 

R TSCA, JS, AND LTJS, RESPECTIVELY . 
ORl<-OR2<-OR3+ORU<-OR5<-OR6<- (NUMSTEPS , (ltplW) )p0 
TR1+TR2+TR3 +TRV+TR5+TR6+ (NUMSTEPS , ( 1 * p VN ) ) p 0 
M 1«-1 

r CALL A5B TO OBTAIN BASE ESTIMATE S3 AND 
R EMPIRICAL ESTIMATE S4 : 

A5B 

r HERE IS THE EXPONENTIAL SMOOTHING FUNCTION 
NEXT : R+ (ALPHA x X4 )+ (S3 x ( 1 -ALPHA ) ) 

R CALL RISKO AND RISKT AND PLACE THE RESULTS 
r FOR THIS ALPHA VALUE IN THE OR AND TR ARRAYS 
RISKO 
RISKT 
M2 + 1 

ONE : OR 1 [Ml : M2 ] +R1 [M2 ] 



OR 2 
OR 3 
OR 4 
OR 5 
OR 6 
TR 1 
TR 2 
TR 3 
TR 4 
TR 5 



Ml 

[Ml 

Ml 

Ml 

[Ml 

"Ml 

Ml 

Ml 

[Ml 

[Ml 



TR6 [Ml 
M2«-M2+l 



M2]*-R2 [M2] 
M2]<-i?3 [M2] 
M2]<-A4[M2] 
M2]«-i?5[M2] 
M2]«-i?6[M2] 
M2 ] +RA0 [M2 ] 
M2] -*rRAT\.M2} 



M2 ] +RM 
M21+RT 
M21+RJ 
M21+RL 



M2] 

M2] 

M2] 

M2] 



->0NE* \ (M2<,2 ) 

NEW STEP : ALPHA+ALPHA +STEP 
M1«-M1 + 1 

+NEXT* \ (ALPHAZUB ) 

V 



Figure A. 1 APL Function ALPHAHAT 

James-Stein shrinkage factors calculated for the base esti- 
mates and the DEE values used in limited translation 
James-Stein for the base estimates are also used for the 
empirical estimates. ALPHAHAT then smooths the empirical 
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Cl] 

[ 2 ] 

[3] 

'4' 

= 5! 
[6: 
[7 3 
[ 8 ] 

[9] 

[10] 
[ 11 ] 
[ 12 ] 
[131 

[14] 

[15] 



V A5B 
DEE+ 0.8 

AN- 3 31 1 +NA 5 
y-e 3 31 1 +YA5 
VN+VNA 5 
7Y<-FyA5 

AW-*- 3 31 1 +AWA5 
SS+(pVN)ol 

ESTIM TO PRODUCE 
1977-1979 DATA : 



>S«- t P 
fl CALL 
R FROM 
ESTIM 
S3+R 
R CALL 
r X4 FOB 
XFOUR 
G+ 0 
V 



BASF ESTIMATE S3 



XFOUR TO 
1980: 



PRODUCE THE EMPIRICAL RATES 



Figure A. 2 APL Function A5B 

estimates onto the base estimates and calls RISKO in Figure 
A. 8 and RISKT in Figure A. 9 to provide the figures of merit. 
RISKO calls BINCONV in Figure A. 10 to produce original scale 
loss rates using the inverse arcsine transformation seen in 
Appendix B for use in the chi-square FOM formula. The 
resulting arrays, which provide figures of merit for all 
estimates and for all validation years, can be analyzed to 
find the a values which minimize transform FOM. 

The above example would find the figures of merit for a 
3-year base period for aviation 1LTS. In order to obtain the 
figures of merit for the other five study cells, functions 
much like A5B are created to set the global variables equal 
to their respective loss and inventory figures. To make 
these calculations for 4- and 5-year base period lengths, 
functions called XFIVE and XSIX are written. These functions 
are the same as XFOUR except that XFIVE finds the empirical 
rates for 1981 and XSIX does so for 1982. The A5B-type func- 
tions are also changed accordingly to alter the number of 
years used to calculate the base estimates. 
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V P+ESTIM; I-,S;U;C1 ; C :K :M:ZB ; ZBBA ;AGO ;AGT ;D 
fl CALCULATES THE AGGREGATE {ORIGINAL ) , 
fl AGGREGATE (TRANSFORMED) , MPE (ORIGINAL) , 
r ISC4 (TRANSFORMED) . JAMES-STEIN , AftP 
fl LIMITED TRANSLATION JAMES-STEIN (LTJS) 
fl ESTIMATORS FOR Z 4S SCREENED BY D. 
fl EACH PAGE OF THE R ARRAY IS A SEPARATE 
fl ESrJMrE. THE DEE WlIt/ES BSEP HEREIN 
fl £BE SEZ 1 1/7 PEE CALLING PROGRAM .A5B , 
fl 48B, ErC.. P JS PEE SCREENING 
fl MATRIX CREATED FROM AVERAGE INVENTORY 
fl WHICH SCREENS OUT THE ZERO INVENTORIES. 
r 



N 

_ ($4ft*0)x0 
>Z))pO 



Z+Y BINPREP 
D*- ( +/ 3 12 
B-<- C 6 , ( 1 ■)' ( 

S+SUNSQ * 

P*(+/,Yx(pY)pP)++/,ftx(pft)pP 

(0 . 5+(+/E* 1 + pft) )*0.5)X lo"l + 2xP 

agt«-(+/,zb)++/,p 

ZBBA+- ( 0 ZB ) p ZEE 

e-e-(CxMxl )+SEJxM=l 

C-e-OT 1-SEJ-e-G 

+(2+l+DLC)xiO<C 

C-«-or l-SEJ-«-(E-3 )* (2-K-K*M)x* /S 

U+ | ( ZB-ZBBA )*SHJ* 0 . 5 



LI 
C 2 
[3 , 
4 : 
5 



^•(l't'pB)p^er 



ft 



«-y WEE 
«-ZB 

<rD*ZBBA+C*ZB-ZBBA 
EE:Cl«-0r 1 - SEE-*- ( 1 LPEE-rE )*SHJ 
B [6 5 ; ]«-PxZBB4+ClxZB-ZBB4 
V 



Figure A. 3 APL Function ESTIM 




\7 Z*<-Y BINPREP N 

n PEEPS PEE FREEMAN-TUKEY VERSION OF THE 
fl 71BC Sift IB4ftS FSB BINOMIAL DATA 
r Y_IS_LOSSESx ft IS INVENTORY 

Ito! 5x ( (0^5+ft)*0. 5 )xZ+ _ lo”l+2x (Y+l )*ft+l 

V 



Figure A. 4 APL Function BINPREP 

Finally, in order to produce the original scale esti- 
mates, the A5B-type functions are modified to perform the 
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AND SSB FOR Z. ALSO 
(ZB) AND GRAND MEAN 
, BOTH DIRECTLY FROM 



V X+SUMSQ Z\ SSE 
[ll n CALCULATES THE SSE 

[2] r CALCULATES THE MLE 
~ 3 j r OR AGGREGATE (ZBB) 
s 4 ] fl TRANSFORMED DATA. 

’ 5 ] K++/ ,D 

= 6 ] ZB+D* (+/Z HM-^l + pZ 

[ 7 ] ZBB<-( + / ,ZB)iK 

[8] X<-+/,P* + /(Z-(pZ)pZSS)*2 , x ^ „ 
= 9] X+(X-SSE),SSE++/ ,Px+/(Z-(pZ)pZB)*2 

V 



Figure A. 5 APL Function SUMSQ 



V Z+Y MLE N-,M 1 

II] R CALCULATES THE MLE IN THE ORIGINAL SCALE 

[2] fl AND TRANSFORMS IT INTO ARCS IN SPACE. 

[3] £><-( + / 3 12 <stAN* 0 )x 0 

[4] Ml*-((p + /Y)pD)x(+/y) + + /jV 

.5] Z«-Dx ( (0 . 5+ (+/N)*1 tpAO*0 ♦ 5 )x lo l+2xMl 



Figure A. 6 APL Function MLE 





V XFOUR 






Cl] 


fl THIS FUNCTION WILL CALCULATE THE TRANSFORMED 


[2] 


n EMPIRICAL ATTRITION RATE X4 


FOR 


USE IN THE 


[2] 


fl EXPONENTIAL SMOOTHING MODEL 


FOR 


ALL SIX 


[3] 


fl ESTIMATORS . CALLS ESTIM. 






[4] 


N+ "1 31 1 +NA5 






[5] 


~1 31 1 +YA5 






C6] 


AN*- '1 31 1 +ANAS 






[7] 


ESTIM 






C8] 


X4<-/? 

V 







Figure A. 7 APL Function XFOUR 

smoothing operation with the a values fixed at their optimal 
levels for each of the six estimators to produce the 
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7 RISKO;D:K;ViBiAR:S;A;NV 
COMPUTES TkE okldNhL SCALE RISK. 

VY AND VN ARE THE VALIDATION 
YEAR LOSSES AND INVENTORY. VN , VY , 

R, AND AN MUST BE IN THE WORKSPACE. 
FIGURES OF MERIT (RISKS) ARE PLACED 
IN VECTORS R 1 THRU R6 . WHICH CORRESPOND 
TO ORIG AGG, ORIG TRANS , MLE , TSCA, 

JS, AND LTjE, RESPECTIVELY . 

'k+ 1 

D+v/AN>0 

V+((pV)oD)x V+BINCONV 
arA(pvy)pd)xvy*vn 

S+(V* O)a(KxI) 

R 1 <rR 2 +R3 +R' 4 *-R 5 <-R 6 «- ( 1 + p VN ) p 0 
LM:A<r(pV)pARlK; :] 

RR^^pF^D^xSxTVVx (Sx (A-V)*2 )*Fx (1-V ) 
i?l IK1<-+/ ,SSZK; ;]x££[l; ;] 

R2 [A'] -*- + /, SSLK-, ; ] *RRL2 ; ; J 
AA1«-+/,SL3 ; ;] 

R3[A]++/,SS[A; 

R4[A]«- + / ,SS ZK ; 

R5ZK3++/ ,SSZK; 

R6ZK1++/ ,SSZK; 

K+K + l 

(Kl + pV«) 

V 



xi?i?[ 3 : 
xRRZ 4 : 
*RR Z 5 
xAR[6 



x(KK++/ ,D)iKK 1 



Figure A. 8 APL Function RISKO 

smoothed transformed figures. Function BINCONV in Figure 
A. 10 is then called, which inverts the transformation and 
yields the original scale attrition rate estimates seen in 
Tables 10 thru 23. 
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Figure A. 9 APL Function RISKT 




V BINCONV:VO ;V1 ;N1 ;D 

A INVERTS ARC SIN TRANSFORMATION . R IS THE 
A ARRAY TO BE TRANSFORMED . N IS THE CENTRAL 
a INVENTORY FOR THE ESTIMATION YEARS . 
D+v/AN> 0 

B-e-0. 5xl + lofi* (0. 5 + (Wl«-(pB)p+/AO*l + pAO*0. 5 
V0+R< (-Ot2)x (M+O . 5 )*0. 5 
Fl«-B>(o*2 )x (Wl+0 . 5 )*0. 5 
B-e-( (pB)pB)xPl+5xF0vFl 
pH 31 p fl 

V 



Figure A. 10 APL Function BINCONV 
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APPENDIX B 

FREEMAN-TUKEY ARCSINE TRANSFORMATION 

1. GENERAL 

Because the TSCA, James-Stein and limited translation 
James-Stein techniques make the assumptions that the distri- 
bution of the number of losses is normally distributed with 
constant variance, and because the binomial model for the 
loss data does not meet these assumptions, a transformation 
is needed. The Freeman-Tukey arcsine transformation produces 
values for which normality and constant variance become more 
tenable assumptions. 

Robinson demonstrated in his thesis [Ref. 2: pp. 74-79] 

that both the normality and variance assumptions are compro- 
mised somewhat for low values of n and p. He therefore 
concluded that the Freeman-Tukey transform is unreliable at 
such values. Therefore, the validity of the results for 
James-Stein estimation and limited translation James-Stein 
must be questioned in this analysis as they were in 
Robinson's [Ref. 2: p. 19] . The following two equations are 

represented in APL by functions BINPREP (transform) and 
BINCONV (inverse transform). See Appendix A. 

2. THE TRANSFORMATION 

The equation for the transformation is: 

x = 0. 5(n+0. 5) 1 / 2 sin _1 ( 2y/( n+1 ) -1 ) ( B. 1) 

+ sin" 1 ( 2( y+l)/( n+l)-l) 

This equation transforms raw losses, y, into transformed 
losses, x using the central inventory, n. 
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3. THE INVERSE TRANSFORMATION 



To invert the transformation and produce the rates in 
original space, use the following set of equations: 



n ij = (1/T)£ n^t), for all i 



(B. 2) 



V 



ij = x i j/( n i j + - 5 ) • 5 



(B. 3) 



0 Vij.le.-nZ2 (B. 4) 

r = . 5( 1+sin Vij(t)) if -n/2<Vi j <n/2 

1 v j_ j . ge. rc/2 

where n^ is the central inventory for the i^ LOS and the 
jth of, i s the corresponding transformed attrition 

figure, and Vij the corresponding loss rate estimates in the 
original scale. 
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APPENDIX C 

ANALYSIS OF OPTIMAL ALPHA OF TRANSFORMED AND ORIGINAL SCALES 
1. GENERAL 

The following tables give the optimal values of a for 
transformed and original scales for the three-year base 
period used in the production of the attrition rate esti- 
mates in Chapter IV. The values of a listed are those which 
produce the minimum figures of merit for validation year 2, 



As one can see from the comparison tables, the a values 
producing the minimum figures of merit for transformed scale 
are very close in most cases to the a values producing the 
minimum figures of merit for original scale, with the 
notable exceptions being MLE for combat support and ground 
combat. Each of the aggregates is discussed below. 



1982. 



TABLE 24 



COMPARISON OF TRANSFORMED AND ORIGINAL ALPHA 
AVIATION AGGREGATE 



1LTS 

AGG ORIG 
AGG TRANS 
MLE 
TSCA 
JS 

LTJS 



TRANSFORMED (X 



ORIGINAL a 



1. 00 
. 00 



1. 00 

1 . 00 

1 . 00 

1 . 00 



1 . 00 
1. 00 
1. 00 
1 . 00 
. 98 
. 98 



LTC0LS 



TRANSFORMED « 



ORIGINAL a 



AGG ORIG 
AGG TRANS 
MLE 
TSCA 
JS 

LTJS 



1. 00 

1 . 00 

1. 00 

1 . 00 

1. 00 

1 . 00 



1. 00 
1. 00 
. 98 
1. 00 
1. 00 
1 . 00 



Note: Optimal a values are for validation year 1982. 
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2. AVIATION 

Looking at Table 24, one sees that for aviation, the a 
values match up almost perfectly except for the transformed 
aggregate for 1LTS, which has values at completely opposite 
ends of the limits of a for the two scales. However, since 
the aggregate transform estimator shows very poor perform- 
ance throughout the analyses in Chapters III and IV, with 
much larger figures of merit than the other estimators, 
inconsistencies like this are not of much concern. 





TABLE 25 








COMPARISON OF TRANSFORMED AND ' 


ORIGINAL ALPHA 




COMBAT 


SUPPORT AGGREGATE 






1LTS 


TRANSFORMED 


a 


ORIGINAL 


a 


AGG ORIG 


1. 00 




1. 00 




AGG TRANS 


. 00 




. 00 




MLE 


. 56 




1. 00 




TSCA 


. 60 




. 50 




JS 


. 90 




. 90 




LTJS 


. 64 




. 56 




LTCOLS 


TRANSFORMED 


a 


ORIGINAL 


a 


AGG ORIG 


. 00 




. 00 




AGG TRANS 


1. 00 




1. 00 




MLE 


. 38 




. 34 




TSCA 


. 34 




. 36 




JS 


. 26 




. 36 




LTJS 


. 30 




. 36 




Note: Optimal a values 


are for validation year 1982. 



3. COMBAT SUPPORT 

From an analysis of Table 25, one sees that the combat 
support aggregate likewise shows a consistency in the values 
of a for transformed and original scales. The only excep- 
tions to this are MLE for 1LTS, which has a difference in a 
.values of a whopping . 44, and, to a much lesser extent, TSCA 
for 1LTS and JS for LTCOLS, which each have a difference in 
« values of . 10. 
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TABLE 26 



COMPARISON OF TRANSFORMED AND ORIGINAL ALPHA 



GROUND 


COMBAT AGGREGATE 




1LTS 


TRANSFORMED 


a 


ORIGINAL 


AGG ORIG 


. 92 




1. 00 


AGG TRANS 


. 00 




. 34 


MLE 


. 94 




. 64 


TSCA 


1. 00 




. 92 


JS 


1. 00 




1. 00 


LTJS 


1. 00 




. 94 


LTCOLS 


TRANSFORMED 


a 


ORIGINAL 


AGG ORIG 


1. 00 




1. 00 


AGG TRANS 


. 00 




. 00 


MLE 


. 42 




1. 00 


TSCA 


. 56 




. 52 


JS 


. 58 




. 50 


LTJS 


. 56 




. 50 



Note: Optimal a values are for validation year 1982. 



4. GROUND COMBAT 

The ground combat aggregate also shows consistency 
between optimal values of a in transformed and original 
scales for all estimators except MLE for both grades and AGG 
TRANS for 1LTS. The lack of consistency between the a 
values for transformed and original scales for MLE seen in 
Table 26 as well as in the 1LTS section of Table 25 is of 
major concern. Unfortunately, using the transformed optimal 
a value to produce figures of merit in the original scale 
for MLE has a big effect on those figures, making them much 
larger. This casts into doubt the ability of the exponen- 
tial smoothing model to produce good maximum likelihood 
estimates of attrition rates. 
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